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ABSTRACT 

We present the first science results from our Hubble Space Telescope Sur- 
vey for Lyman limit absorption systems (LLS) using the low dispersion spectro- 
scopic modes of the Advanced Camera for Surveys and the Wide Field Camera 
3. Through an analysis of 71 quasars, we determine the incidence frequency of 
LLS per unit redshift and per unit path length, i{z) and i{X) respectively, over 
the redshift range 1 < z < 2.6, and find a weighted mean of i{X)= 0.29±0.05 for 
2.0 < z < 2.5 through a joint analysis of our sample and that of Ribaudo et al. 
(2011). Through stacked spectrum analysis, we determine a median (mean) value 
of the mean free path to ionizing radiation at z = 2.4 of X^J^l = 243(252)/i72^ Mpc, 
with an error on t he mean value o f =b43/ tyr,^ Mpc. We also re-evaluate the esti- 



mates of A^£p from lProchaska et al.l (120091 ) and place constraints on the evolution 
of A^£p with redshift, including an estimate of the "breakthrough" redshift of 
z = 1.6 

Consistent with results at higher z, we find that a significant fraction of 
the opacity for absorption of ionizing photons comes from systems with A^hi< 
]^gi7.5pj^-2 ^j^j^ value for the total Lyman opacity of t^^™^"^ = 0.40 ± 0.15. 
Finally, we determine that at minimum, a 5-parameter (4 power-law) model is 
needed to describe the column density distribution function /(A^hi, X) at 2; ~ 2.4, 
find that f{Ni{i,X) undergoes no significant change in shape between z ~ 2.4 
and z ~ 3.7, and provide our best fit model for f{Niii,X). 
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Introduction 



The importance of Lyman limit systems (LLS), those quasar absorption hne systems 
with A^Hi > 10^'^'^ cm~^ has been well understood. Studies of LLS at redshifts z > 2.6 
stretch back many dec ades, with th e LLS being some of the first absorption line systems 
studied quantitatively iTytlerl (119821 ) . Through to the present day, ground based s urvey s 
of the LLS span the full ran g e of observab i hty, 2 . 6 > z > 6 (e.g. Sargent et al.l (119891). 



Storrie-Lombardi et all (119941 ). IPeroux et al.l ( ]200l[ ). ISongaila fc Cowid (120101). 



Prochaska et al. 



(2010)). The 1990s also introduced studies of the LLS from space (e.g. IStengler-Larrea et al. 
( 1l995l )). with the largest datasets being provided by the the Faint Object Spectrograph and 
Space Telescope Imaging Spectrograph onboard the Hubble Space Telescope spanning the 
range 0.3 < z < 2.6 (see 
peak in sensitivity a.t z ^ 



Ribaudo et al.l ( 120111 ) for a summary). The space-based studies 
1, declining to both higher and lower redshifts. Nearly all early 
studies of the LLS provided a general picture of a rapidly evolving population with an inci- 
dence frequency described by a (1 + z^ power law, and 7 ~ 1.5. Early studies also provided 



e.g. 


Tvtler 


( 


1987 


); 


Lanzetta 


( 


1991 


); 


Peroux et al 


J( 


2003 


); 


O'Meara et al. 



( I2OO7I ): iProchaska et al.l (I2OIOI )) have considered more complicated descriptions. 



The optical depths from LLS at 1 Rydberg, combined with their increased frequency 
per unit redshift compared to the higher A^hi damped Lyman alpha systems (DLAs; A^hi 
> 10^°'^cm~^), suggest that they strongly influence, if not dominate the attenuation of 
photons emitted by galaxies and quasars, and thus set the intensity of the extragalactic 
UV background, and determine the mean fre e path of ionizing photons in the intergalac- 



tic medium (IGM; e.g . jMiralda-Escude (120031 ) : IProchaska et a. 



DalPAdio et al.l ( 120081 ) 



Calverley et al. 



(I2OO9I), hereafter PWO 09; 



j201ll ) JMcQuinn et al.N201l[ ): Iflaardt fc Madaul ( l2012[ )) 



Thus, LLS play a crucial role in the cosmological history of interplay between radiation and 
baryons. 

Both the low density, highly ionized Lyman-a forest ( A^hi < lO^^'^cm^^ ) and the high 
density, predominantly neutral DLAs have been well studied and constrained in part because 
of their relative ease of observation. The Lyman-a forest is ubiquitous in quasar spectra, 
is visible from the ground at z > 1.6, and is relatively easy to constrain in terms of A^hi 
owing to its lack of saturation in the Lyman series. One must go to higher order Lyman 
series transitions to determine t he A^ht as A^ht incr eases, but this can be achieved with high 
S/N, high- resolution data (e.g. iRudie et al.ll2012l ). The DLA show strong damping wing 
features in their Lya line, al lowing for easy determination of their A^hi even in low-resolution 
data (e.g. IWolfe et al.ll2005l ). By contrast, the LLS, although frequently observed, are by 
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comparison poorly constrained, as both the Lya line and the Lyman break must be well 
covered in the spectrum to constrain the A^hi (e.g. iProchter et al.ll2010l ). At the low end of 



the LLS H I column density range, significant wavelength coverage at rest wavelengths in 
the LLS frame of A < 912A are required to accurately determine the A^hi through the shape 
of the recovery in fiux. These complications place stronger constraints on the redshift range 
where large statistical surveys for LLS can be done from the ground, namely z > 2.5. 

Nevertheless, our motivation for studying the LLS is clear. A complete census of ab- 
sorption across the full range of A^hi is required to understand the properties of the IGM 
at a given redshift, and the evolution in redshift of that census provides tests of models of 
structure formation, evolution i n the UV background, a nd the interplay between galaxies and 
their environment. At z ~ 3.7. iProchaska et al.l ( 120101 ) (hereafter POWIO) analyzed spectra 
from the SDSS DR7 catalogue to establish the incidence frequency of LLS and to characterize 
a number of biases inherent in LLS analysis at any redshift. PWO09 also use the SDSS DR7 
catalogue, and produce a measurement of the mean free path to ionizing radiation, A^£p, 
over the same redshift range in POWIO. Combined with the incidence frequency measure- 
ments and measurements of the Lyman-a forest, LLS, and DLA in the literature, POWIO 
determined the column density distribution function f{Niii,X) at z = 3.7. This /(A''hi,X) 
underscores the importance of the LLS on our complete understanding of absorbers in the 
universe, as it showed that a number of strong infiections in f{Ni{i,X) exist at this red- 
shift. Recent dev elopments in theoretical models and simulations of the LLS regime (e.g. 
Altay et al.ll201ll ) have made significant improvements in reproducing these infiections with 
the inclusion of self shielding. The path forward is now clear: to move beyond a simple 
counting of the LLS and to extend the knowledge of /(A^hi, X) and A^£p to both higher and 
lower redshifts to provide a determination their evolution with cosmic time, and to better 
understand their physical nature. 

At lower redshifts {z < 2.5), a determination of f(Nm,X) and A^fp which includes 
the contribution from LLS is made difficult by the simple fact that one must go to space 
based telescopes to obtain coverage of the Lyman break at rest wavelength A = 9 12 A. The 
region 1 < z < 2.5 is particularly difficult, as most space-based instruments have poor NUV 
throughput at high resolution. Previous studies of LLS at 2; < 2.5 offered d i ffering results 
regarding the evolution of the incidence frequency, i{z) of the Us. iLanzetta dim) argued 
for a fairly constant £(z) with rap i d evo lution at higher z, whereas IStengler-Larrea et al. 



(119951 ) and IStorrie-Lombardi et al.l (Il994f) argued for a single power law evolution over the 
entire range < z < 4. Recently, iRibaudo et al.l ( 1201 ll ) performed an HST archival study 
of the LLS frequency distribution at z < 2.5, with the bulk of their statistical power at 
0.75 < z < 1.5. They find that the LLS frequency with redshift i{z) is well described by 
a power-law i{z)(x (1 + z)"^ with 7 = 1.33 ± 0.61 at z < 2.6. They place constraints on 
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f{Ni{i,X) and A^fp (finding that the former is incompatible with a single power-law), but 
do so largely by relying on measurements at different redshifts than the bulk of their sample 
for the higher column density LLS and the DLA. 

In this paper, we present the first set of scientific results from our campaign to better 
explore the LLS at redshifts where our knowledge of optically thick systems is poorly con- 
strained, namely 1 < z < 2.6. Our data sample comes from a survey performed with the 
Hubble Space Telescope using the low-resolution NUV spectroscopic modes of the Advanced 
Camera for Surveys and the Wide Field Camera 3 instruments. The survey is described in 



the first paper in this series lO'Meara et al.l ( 120111 ). which we refer to hereafter as Paper I. 



This paper is organized as follows: Section |2] summarizes the results of Paper I, Section [3] 
describes how we model LLS absorption in our sample. Section H] describes our statistical 
analysis and results for the incidence frequency of the LLS, Section [S] describes our analysis 
of stacked spectra and the determination of the mean free path of ionizing photons under 
specific assumptions regarding the stacked QSO SED and any intrinsic infiections it might 
have. Section [6] presents our analysis of the column density distribution function, and Section 
[7] offers a summary and discussion. 

Unless stated otherwise, all results in this paper assume a "standard" Lambda-|-CDM 
cosmology with Hq = 72km s~^ Mpc~^, Qm = 0.26, and = 0.74. 



2. Summary of Quasar Selection and HST Datasets 



In Paper I, we reported in detail the quasar sample and i/ST observations for our survey. 
We also described the reduction procedures and presented the extracted ID spectra that form 
the basis of the following analysis. This section provides a brief summary of Paper I. 

For our i/S'Tsnapshot programs in Cycles 15 and 17, we selected 100 quasars at Zem ~ 2.5 
from the SDSS Data Release 5 based primarily on their optical photometry. Specifically, we 
restricted the list to quasars with g < 18.5 mag, 2.30 < z^m < 2.60, and spectra without 
very strong associated absorption. B ecause the SDSS qua sar spectroscopic sample is based 
primarily on optical color-selection (IRichards et al.l 120021) . it is possibl e it may be biased 
relative to a complete quasar sample (e.g. IWorseck fc Prochaskal l201ll ). Indeed, PWO09 
demonstrated that the cohort of SDSS quasars at 2; ~ 3.5 are biased toward sightlines with 
strong Lyman limit absorption in the M-band. For our selection criteria, the presence of an 
intervening LLS should have negligible effect on the quasar colors. Nevertheless, quasars at 
z ~ 2.5 have a typical color that is sufficiently close to t he stellar locus that the SDSS team 
chose selection criteria which favor UV-excess quasars. IWorseck fc Prochaskal (120111 ) have 
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demonstrated that this biases the spectroscopic sample at z ~ 2.6 to have bluer spectral 
energy distributions (SEDs) than a complete sample (see their Figure 16). We return to this 
issue later in the paper. 

From this list (see Table 1 of Paper I), 18 quasars were queued in Cycle 15 to observe 
with the ACS/PR200L prism and 53 quasars with the WFC3/UVIS-G280 grism in Cycle 17. 
We developed customized software to extract, wavelength calibrate, and flux the data to 
produce a fully calibrated ID spectrum for each quasar. The ACS/PR200L spectra cover 
observed wavelengths A = 1500 — 5000 A and the WFC3/UVIS-G280 spectra span roughly 
A = 2000 — 6000 A. Our spectra have relatively high S/N per pixel down to A ~ 1800 A 
for the ACS dataset and A ^ 2000 A for the WFC3/UVIS-G280 dataset. This reflects the 
relatively high UV fluxes of the quasars and the low dispersion of the spectrometers: at 
A = 2500A the spectra have FWHM 60A. Uncertainties in the wavelength calibration are 
on the order of 2 pixels, in the form of a rigid shift in pixel space. Table [1] summarizes the 
quasars comprising the survey. 



3. Fitting for Lyman Limit Absorption 

The principle goal of our i!f5T program was to survey z ~ 2.5 quasars for Lyman limit 
absorption at z < 2. At these redshifts, the Lyman limit lies below Earth's atmospheric cutoff 
and one requires space-borne, UV spectroscopy. Because the Lyman limit is a continuum 
opacity, its analysis does not require high spectral resolution. Instead, one prefers well-fiuxed 
spectra with high S/N, and, ideally, quasars without complex continua. Our dataset nicely 
achieves these criteria (Paper I). 

Our approach to the Lyman limit analysis is straightforward. First, we must estimate the 
observed flux of each source in the absence of Lyman continuum opacity. This quantity is not 
the intrinsic flux of the quasar; the Lya forest scatters light at rest wavelength A^ < 1215A 
beginning with Lya opacity and eventually including the full Lyman series. Because of the 
low spectral resolution of these spectra, we cannot identify individual absorption lines and 
therefore observe the intrinsic quasar flux attenuated by a relatively smooth Lyman series 
opacity (see § 15.11 for a full description). Therefore, our approach was to flt a template 
spectrum to the flux at Ar ~ 9 50 — 1150A. We have found that the data are reasonably well 



described by the iTelfer et al.l (120021 ) radio-quiet quasar template spectrum fj^^'^^^'^ shifted to 



the observer frame and allowing for a scaled normalization C and power-law tilt a, i.e., 

rrw < 1200A) = c f^- ( ^ y . (1) 



- 6 - 



We emphasize that even though the Telfer spectrum is intended to represent the mean 
intrinsic SED of 2; ~ 1 quasars, fx°^^^ represents the intrinsic flux of each of our quasars 
attenuated by Lyman series opacity. The values for C and a were determined using a 
custom, interactive GUI that allows one to visually compare a model continuum with each 
quasar spectrum. In general, the data at Arest = 950 — 1150A in the QSO rest-frame constrain 
the C and a parameters with reasonable confidence (~ 10 — 20% uncertainty). The analysis 
on each quasar was done independently by JMO and JXP, and we report our preferred 
models in Table [H We then extrapolate this model to Ar < 950A. 

With Z™"^*' estimated as above, we proceeded to fit any significant and sudden drop in 
the flux below Arest = 912 A as Lyman limit opacity. Specifically, we model the observed 
flux as the continuum flux described above modulated by one or more "systems" yielding 
appreciable Lyman limit opacity. Each system is characterized by a redshift Zabs and an 
optical depth at the Lyman limit Tg^ ~ cm~^. In the following, we refer to 

systems with Tg^ > 2 as Lyman limit systems (LLSs) and systems with Tg^ < 2 as partial 
Lyman limit systems (pLLsfl), although we occasionally use LLS to refer to all systems 
exhibiting detectable Lyman limit opacity. The resulting model is 



Y-modcl rconti 

Jx — J X 



exp 



N 



LL.i 



(2) 



where the sum is over all absorbers identi fied along the sigh tline (zabs < -2^em) and tx is the 
opacity in the rest frame of the absorber (jVerner et al.lll996l ). 



In a minority of cases (~ 30%), the observed spectra are well described by one or zero 
Lyman limit systems. For the remainder of cases, however, the data are best, and very 
well, modeled with multiple LLS with Tg^ varying from 0.2 to 2. As with the continuum 
estimation, the modeling was performed by hand using a custom GUI that accounts for 
the spectral resolution of the HST data. We also experimented with x^-minimization algo- 
rithms, but these gave unrealistically small statistical errors («i 1 — 2%) owing to the high 
S/N of the spectra (and high reduced values because the models did not include Lya 
forest absorption). Indeed, the results are dominated by the systematic error of continuum 
placement and, to a lesser extent, line-blending by the IGM. For these reasons, we have 
proceed with a by-hand approach and estimate the uncertainties through inter-comparisons 
amongst multiple authors, and further with the results from a Monte Carlo analysis. 



^We note the difference between pLLS as defined here, and PLLS, which is a "proximate" LLS, i.e. one 
within 3000 km s^^ of the QSO redshift. 
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Each spectrum was modeled independently by JMO and JXP and then the authors 
compared these results and agreed, along with a third author, GW, on a solution for cases 
in dispute. For the systems with > 2 listed in Table [31 the two authors had agreed in 
35/39 cases there was a rc[f^ > 2 system within 6z = 0.05 of the redshift listed in the table. 
Of the 4 disparate cases, most had Tg^ ^ 2. In these cases, the redshifts have an RMS of 
0.01. For the systems in the table with 0.5 < < 2, both authors agreed on a system 
within 6z = 0.05 in 20/27 cases. The RMS in the optical depths between the two sets of 
models is ~ 0.05 and the RMS in redshifts is ~ 0.02. 

To gauge our ability to successfully identify LLS absorption, we performed an analysis 
on 100 mock quasar sightlines. The sightlines were generated by drawing randomly from 
Monte Carlo absorption line lists, each created using the best-fit column density distribution 
function described in 16.41 Absorption lines were generated from these lists, convolved to the 
WFC3 resolution, and were placed on absorption-free QSO spectra with varying power law 
tilts and signal to noise ratios chosen to best mimic the data variety of our sample. The 
spectra were then surveyed by one author (JMO) for LLS absorption using the same interface 
as with the real data. The results of this analysis are as follows: When we impose the survey 
definition criteria given in l4.H the mock LLS spectra contain a total of 91 Tg^ > 0.5 absorbers 
in the survey path, with 72 systems having Tg^ > 2, 11 having 1 < Tg^ < 2, and 8 with 
0.5 < 

Tgi2 — l-O- JMO identified 76 systems having > 2, 9 having 1 < Tg}^ < 2, and 
10 with 0.5 < Tg^i^a < l-O- The primary nature of the discrepancy between the input mock 
LLS sample and those recovered dealt with systems very near to the optical depth boundary 
of each bin. For example, three of the four Tg^ > 2 discrepant cases, the input mock 
spectra had 17.4 <logA^Hi< 17.5, but were given a column density of either logA''Hi= 17.5 
or logiVHi= 17.55 Likewise, one of the two 1 < Tg^ < 2 discrepancies and two of the three 
0.5 < ^912 < 1-0 cases stemmed from the mock LLS lying within logiVni— 0.05 of the 
boundary between bins in Tg^ . In one case, a single Tg^ > 2 system was used to describe 
two mock LLS with Tg^ > 1 which lied within Sz = 0.1 of each other. In only a single case 
was a completely incorrect system used to model the input mock spectrum, where a spurious 
Tg^ ~ 1 system was included. Although not included in any of our statistical analyses below, 
> 50% of all 0.2 < Tg^ < 0.5 absorption was correctly identified, confirming our ability to 
be sensitive to such absorption. 

Figure [1] gives a further exploration of the range of uncertainties in our data. In the 
upper panel, we show the estimates of log A^hi in LLS within the range log Ahi> 16.9 derived 
from fits to the mock spectra compared to their input values. We choose this range to fully 
explore the range of values for which we can make more precise statements about A^hi, 
namely the statistical survey range 0.5 < Tg^ < 2.0 with additional measurements at higher 
optical depths when the data allows for an A^hi do be determined in the cases where we 
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can observe the recovery in flux below the Lyman limit (although we note that all > 2 
absoprtion is lumped into a single bin for our statistical analysis). We find overall good 
agreement between the fit values and the input mock values for A^hi over this range with a 
mean difference between fit and input of —0.02, and a = 0.12. In the lower panel of Figure 
[H we explore our ability to properly determine the redshift of the LLS for the Tq^^ > 0.5 
LLS. We find a mean difference between fit and input 2;iis of —0.002 and a = 0.014, again a 
good agreement. 

These comparisons provide an estimate of the systematic uncertainty related to the LLS 
modeling. In addition, we estimate the systematic uncertainties of continuum uncertainty 
and system blending to be 1% f~ 1000km s~^) for the redshift estimation and ~ 15% for 
the optical depth measurementso Figures |2] and E] present the spectra and adopted models 
for WFC3/UVIS-G280 and ACS/PR200L respectively, and Table tabulates the results. 



4. The Statistics of LLS at z ^ 2 

4.1. Survey Definition 

From the spectral analysis of the previous section, it is relatively straightforward to 
perform a statistical analysis of Lyman limit absorption. A fundamental description of 
IGM absorption systems is the iVni frequency distribution, /(A^hi,-^)- Akin to a luminosity 
function, /(A^hi? X) gives the number of absorbers with column density A^hi per ciA^ni interval 
and per path length dX, where 

dX=^{l + zfdz . (3) 
H [z) 

The absorption path length was introduced to yield a constant f{Niii,X) for systems that 
arise from discrete absor bers (e.g. galaxies) with a constant product of comoving number 



density and physical size ( IBahcall fc Peebleslll969l ). Standard practice is to survey quasars for 



absorbers as a function of H I column density and estimate f[N}ii,X) from simple counting 
statistics. 

Owing to uncertainties in the continuum placement and limitations in S/N limit, analysis 
of our quasar spectra has a limiting sensitivity to Lyman limit absorption of Tg^ > 0.2 or 
A^Hi ^ lO^^'^cm"^. However, as explained below, we only include those systems with Tg^ > 
0.5 in our statistical exploration of incidence frequency and column density distribution. 



^AU models assume b = 25km s ^; this has no bearing on our resuUs given the low spectral resolution. 
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Furthermore, it is impossible to yield any estimate of Tgf^ beyond a lower limit for systems 
with > 4 because all such systems have zero measurable flux below the Lyman limit. 
Although we can use the recovery in flux below the Lyman limit to obtain Nm for LLS with 
2 < r^^2 < 4, such cases require significant spectral coverage which is not always available in 
our data. As a result, our statstical sample will consder all systems with > 2 together 
as a single population. § In short, the analysis has sensitivity to Tg^ over approximately 
one order of magnitude and many systems only yield a lower limit to Tq^. Therefore, one 
may estimate f{Niii,X) directly in only a narrow range of H I column densities. For these 
reasons, studies of Lyman limit absorption have provided rather limited direct constraints 
on f{Nni,X). 

Instead, researchers have focused primarily on the zeroth moment of /(iVHi,X), i.e. the 
incidence of Lyman limit absorption per unit path length 



^(^).>rn..= j f{Nni,X)dNni (4) 

-'VhI, limit 

or per unit redshif10 i{z). The incidence of Lyman limit systems is measured from a limiting 
optical depth rnmit corresponding to a limiting H I column density A^ni.iimit, generally set by 
the spectral quality. 

In principle, the i{z) quantity is a direct observable. Standard practice is to estimate it 
from the ratio of the number of systems detected within a redshift interval A/lls to the total 
redshift path surveyed Az: 



LLS 



(5) 



Observation ally, the r e are several factors that complicate such an LLS survey (see e.g. Ap- 
pendix C in iTytlerl ( 119871 )). First, it is very challenging to successfully identify multiple 
LLSs that have small separations in redshift space. One is limited by the spectral resolution, 
the S/N of the data, the precise characteristics of the systems, and the density of the Lya 
forest at the wavelengths of interest. For our data, we cannot resolve two LLS located within 
\6v\ < 10,000 km s~^ {6z ^ 0.1) and therefore consider all such complexes as a single LLS 
with a summed optical depth. Second, the presence of a partial LLS reduces the quasar flux 



■^One can improve the estimate for iVni from the Lyman series, e.g. the damping wings of Lya, but our 
spectral resolution precludes such analysis. 

''Also commonly referred to as dAf/dX and dAf/dz. 
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and the resultant spectral S/N but may not prohibit the search for an additional LLS. As 
described in detail in POWIO, one must establish a strict and proper criterion for the survey 
path or the results will suffer from a "pLLS bias" . 

To this end, we establish these criteria for the survey: 



1. The ending redshift Zg^id is set to be 3000km s~^ blueward of the quasar emission 
redshift. This is established in part to account for the relatively large uncertainty in 
quasar redshifts and also to minimize any biases from gas in the environment including 
the quasar. 

2. If there are no LLS detected, the starting redshift ^start = (2200 A/911.7641 A)-l = 1.4 
for the WFC3 data and Zstart = 1-2 for the ACS data. This provides enough spectral 
coverage (~ 200 A) to confidently identify an LLS system before the spectra deteriorate 
in quality. 

3. For the WFC3 sample, identification of an LLS with > 0.5 and z = Zabs precludes 
the search for any LLS with z < z^hs and Tg^ < 1. 

4. For he WFC3 sample, identification of an LLS with > 1 and z = z^hs precludes 
the search for any LLS with z < z^hs and Tg^ < 2 . 

5. In both samples, identification of an LLS with Tg^ > 2 and z = z^hs precludes the search 
for any other LLS with z < Zabs-Because of potential sky subtraction complications 
(see Paper I), we have only performed a survey for Tg^ > 2 LLS in the ACS/PR200L 
spectra. 



Together, the final four criteria establish -Zstart values for a range of limiting optical depth. 

With these survey criteria, it is straightforward to evaluate the survey path of our 
dataset. This is described in Figure HJ which presents the so-called sensitivity function g{z) 
for the survey, i.e. the number of unique quasar spectra at a given redshift where one can 
search for an LLS to a given optical depth limit. The survey path is similarly summarized in 
Tables |3] and m where we hst the start and ending redshift for the LLS survey of each Q SO. 



We note that the sensitivity function is comparable to that from iRibaudo et al.l (l201ll ) at 
the low {z ~ 1.5) redshift end of our sample, increasing to nearly twice their value at higher 
redshifts, making the two samples highly complementary. None of the quasars that they 
surveyed are in our sample. 
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4.1.1. Results 



One may now simply evaluate Equation |5] to derive estimates for the incidence of LLS 
absorption at z < 2, as a function of limiting optical depth. Figure [5] presents the results for 
the WFC3/UVIS-G280 survey for rg^^ > g.s, i and 2, and for the ACS/PR200L survey for 
^912 > 2 LLS. We present the results for two redshift intervals, z = [1.2, 2.0) and [2.0, 2.6). As 
one predicts for an expanding universe, we observe a decreasing incidenc e between the two 
redshi ft bins for all samples. Overlayed on Figure [5] are the results from ISongaila fc Cowie 



mm for 

T912 > 1 deter mined over the redsh ift range < z < 6. Our results are consistent 



with theirs and those of iRibaudo et al.l (120111 ). but we note that our survey provides better 
sensitivity for the epochs 1.5 < z < 2.5. 



Equation |3l along with the transformation 

e{X)dX = i{z)dz 

and the cosmological evolution of the Hubble constant 



H{z) = Ho[nA + nrn{i + zy] 



3 1 1/2 



(6) 



(7) 



allows for the determination of i{X) from our sample. We use the same bins in redshift 
and sub-sampling in as for £{z). The results for i{X) are sho wn in Figure El Also 
shown in Figure El is the single power law fit for > 2 LLS from iRibaudo et al.l (1201 if ) 
which combines their analysis of 2; < 2.5 HST spec tra and the z > 3.5 S PSS results from 
POWIO. Our results generally agree with those of IRibaudo et al.l ( 1201 ll ). namely that at 
redshifts z < 2, £{X) evolves only very weakly (in fact, the WFC3/UVIS-G280 results are 
more consistent with no evolution), transitioning to a sharper evolution with z aX z > 2. We 
further concur with their observation that the 2.5 < z < 3.5 span must be fully explored 
to best understand this transition in evolution of i{X). Additionally, Ribaudo et al. (2011) 
use their results to constrain the differential column density distribution, /(A^hi?-^)- We 
address this quantity below, using the additional information gained from an analysis of 
stacked spectra. 

We now pl ace our results for £{ X) from WFC3 /UVIS-G280 alongside th ose obtained 
from the SLLS Jo'Meara et al.lboOTh and the DLA JProchaska fc Wolfel [2009[ ) in Figure H 
Quantitatively, we find that the SLLSs and DLAs (i.e. systems with A^hi > lO^^cm"^) 
contribute a significant fraction of £(X)^>2. By inference, this implies that there are relatively 
few LLS with A'hi = 10^^-^ - 10^^ cm'^ and therefore that /(A^h^AT) is shallow at these 
column densities. We return to this issue in § 16. 1[ We also note with interest the relative 
contributions to (!-{X) that each bracket in r^^^ makes. Specifically, we see that the integral 
contribution to £{X) of the LLS with 0.5 < Tgi2 < 2 is nearly equal in magnitude to that 
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for the LLS with > 2. As i{X) corresponds to the integral of f{Ni{i,X) over iVni, this 
imphes a pronounced change in shape in f{Niii,X) for < 2. We explore this result 
in greater detail in § 16.21 Finally, we combine our results for > 2 for ACS/PR 200L, 
WFC3/UVIS-G280, and the results from the archival study of lRibaudo et al.l toilh . The 
weighted mean for these results at z ~ [2, 2.5] is i{X)r>2 = 0.29 ± 0.05. We adopt this 
measurement in all analysis that follows. Table |5] summarizes all of these results. 



4.2. Systematic Errors in the Incidence Frequency 



To assess the effects of uncertainty in the A^hi and ^abs values for the LLS, we performed 
a Monte Carlo analysis where we randomly modified each system with a normal deviate 
and assuming aN^^ = 0.05 dex and = 0.02 (see § [3]). We then recalculated i{z) 1000 
times for each redshift bin. We find that none of the trials have an i{z) value that lies 
outside the 85% confidence level given by Poisson uncertainty. Furthermore, we explored the 
effects of proximity to the quasar by varying our upper wavelength cut. This simultaneously 
addresses issues associated either with inaccurate quasar emission redshifts, Zqso, or with the 
possibility that quasars arise in overdense regions, which might effect our statistics. To do 
so, we expanded the proximity region from 3000 km s~^ blueward of the quasar emission 
redshift to 10000 km s~^, and find our results to be essentially unchanged. 



In iProchaska et al.l ( 120101 ) we identified and described at length a systematic bias to 
LLS surveys associated with pLLS. This bias occurs when one truncates the search path 
for an LLS owing to the presence of one or more pLLS, as we have done in this paper for 
systems with < 2. It is possible, therefore, that we have over-estimated the incidence of 
systems with Tg^ ~ 1. The pLLS bias, however, is proportional to the incidence of Lyman 
limit absorption which at these redshifts is small. Indeed, we estimate that this is a < 10% 
effect and therefore within the statistical error associated with our small sample size. 



The other, and unavoidable, bias described by IProchaska et al.l ( 120101 ) is related to 
the blending of Lyman limit systems. One or more systems localized to 5z ~ 0.1 will be 
unresolved by our spectral resolution and therefore treated as a single LLS. This raises the 
incidence of Tg^ > 2 systems at the expense of pLLS. Again, this bias is proportional to the 
overall incidence of LL S and at these redshifts the effects will be small (< 10%; Figure 7 of 
Prochaska et al.lboioh . 



- 13 - 



5. Stacked Spectrum Analysis 

In this section we generate an average quasar spectrum from tlie WFC3/UVIS-G280 
spectra and then analyze this dataset to estimate the mean free path to ionizing radiation 
A^fp, following PWO09. We begin with a review and expansion of their formalism. 



5.1. Formalism 

In PWO09 we introduced a new technique for constraining the mean free path to ionizing 
radiation, A^fp. The approach analyzes a 'stacked' quasar spectrum, constructed from a strict 
average of a cohort of quasars with common emission redshift, at wavelengths blueward of 
the Lyman limit. This stacked spectrum, by definition, yields the average observed flux 
yobs q£ quasars at a given redshift. For wavelengths redward of Lya, the flux is dominated 
by the spectral energy distribution (SED) of quasars, with only a minor attenuation from 
the average metal-line absorption of gas in the IGM, t^^^'^^^{z). Blueward of Lya, e.g. 
at ~ II25A, the observed flux fii2^ is the product of the average intrinsic quasar SED 
fiui the average absorption from the Lya forestlfl exp(— r^°) at the redshift ZLya = 
1125A (1 + Zem)/1215.67A - 1, 



/ll25(^em) — /n25 ^^P ''cff (^Lya) 



-Lya/ 



As one progresses to shorter wavelengths (l ower z^^^), T p^g" is observed to decrease because 



the incidence of Lya lines decreases (e.g. iRauchI Il998h . However, one eventually incurs 



additional opacity from Ly/3 from gas at redshift zi^yp = A/1025.72A — 1, and then Ly7 
opacity, and eventually the full Lyman series. 



00 



ri,— (ze„.,A) = 5^r:s(zO , (9) 



n=2 



with Zn = A(l + Z em)/Xn — 1 and n = 2,3,4,... corresponding to Lya, Ly/3, Ly7, etc (e.g. 



Madau et al. 



19961). 



At Ar < 9 12 A, the emitted photons also experience the continuum opacity of Lyman 



^Observationally, it is common practice to assess T^g" to a limiting H I column density, e.g. to exclude 
DLAs or even LLSs. 
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limit absorption r^g {zem, Zgu) with 



2^912 



A, 



, 1 ~l~ Zcxn I 1 



(10) 



912 



This quantity is related to the mean free path A^fp, defined to be the distance that a packet 
of photons travels before suffering an e~^ attenuation (on average). Specifically, we define 
-^mfp to t)e the physical distance from Zem to the redshift z^f* where T^(2;^f', Zem) = 1- Note 
that with this definition, A^^, applies to photons with energies greater than 1 Ryd at 2;em- 
The mean free path to i onizing photons establishes the rnean i ntensity of the extragalactic 
UV background (EUVB; iHaardt &: Madaulll996l : iMeiksinl l2009l ) and, in essence, defines the 
epoch of Hydrogen reionization (i.e. A^fp has negligible value in a neutral universe). 

Estimations of (and thereby A^jp) have been made previously by integrating eval- 
uations of the H I frequency distribution /(./Vht, z) of absorption systems in the IGM (e.g. 
Madau et al1ll999l : IPaucher-Giguere et al.lSoOSah : 



rl^{z9i2, Zera) = j j /(iVffl, z'){l - exp [ -iVHiaph(^') ]]dNmdz' (11) 

»)12 

where Uph is the photoionization cross-section evaluated at the photon frequency. This 
approach is fraught with great uncertainty because: (i) the frequency of absorbers with 
roughly unit optical depth at the Lyman limit t^^2 ^ 1 is very difficult to ascertain; and 
(ii) the incidence of optically thick absorbers (r^]'^ > 2) suffers from significant systematic 
uncertainty (POWIO). 

In PWO09 we defined a new formalism to evaluate r^^^ from stacked quasar spectra. 
We introduced an effective opacity from Lyman limits kll which evolves with redshift and 
frequency as 

/ jy \ -2-75 / ]^ _^ 2 \ -2.75 

kll{z,i^) = HgMz) I I = ^912(2;) I — I , (12) 

V ^^912 / V 1 + ^912 / 

where the second term accounts for the photoionization cross-sectior0 (accurate to within 
1% to Ar ~ 600 A). One may relate this opacity to the optical depth through the standard 



^In PWO09 we assumed a, 1/ ^ dependence which is not as accurate. Furthermore, there is essentiaUy no 
frequency dependence on the effective opacity for systems with very high A^hi, which our approach ignores. 
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definitioiill 

rLL(r,z/)= / «:LL(r',z/)rfr' , (13) 







and relate distance to redshift via cosmology 

dr c c/Hq 



, = (U) 

dz ii + z)H{z) (i + z),/njTTWTIhi 

In PWO09 we performed this analysis at 2; ~ 4 where the Hubble parameter H{z) is domi- 
nated by matter and one can approximate dr/dz oc (1 + z)~^^'^. At the redshifts of interest 
here (2; < 2), we find that a power-law representation, 



dr , ^ dr , 

-(.<2.6)«-(. = 2,6) 



l+Z 

3.6 



-2.315 



(15) 



is accurate to within 1% for z = 1 — 2.5. We adopt this power-law approximation in the 
following, taking dr/dz{z = 2.6) = 322/172^ Mpc. Altogether our expression becomes 

re'^F^(^912,^cm)= 6237Mpc(l + Z912)'-''y" ^^912 (^') ( 1 + C?^' ■ (16) 



^912 



Note that this formalism is primarily introduced to parametrically solve for X^^^ by 
establishing the redshift where t^q = 1. In contrast to the PWO09 analysis of SDSS spectra, 
the WFC3/UVIS-G280 data have the wavelength coverage and sensitivity to produce a 
stacked spectrum to ^ 6OOA. In part, this reflects the larger mean free path to ionizing 
radiation at z ~ 2 than z ~ 4. If one wishes to analyze the stacked spectrum to these 
wavelengths, however, it is not safe to assume that the modulations in flux at Ar < 912A 
are due solely to Lyman limit opacity (as was done by PWOOslfl). First, one must consider 
evolution in the intrinsic SED ySED ^y^^ the s e wav elengths. For example, the 2; ~ 1 average 



quasar spectrum produced by iTelfer et al.l ( 12002| ) follows a power-law /|^° oc A at 



Ar < I200A which implies a 20% higher flux at Ar = 6OOA than at the Lyman limit. Second, 
the Lyman series opacity T^g'^^^ is certain to decrease with redshift. Therefore, the observed 
flux should increase at Ar < 900A in the absence of Lyman limit opacity. 



^This formalism wrongly ignores stimulated emission, but that has a negligible effect on our model. 

^This assumption was justified, in part, because the analysis was performed over a very short wavelength 
(redshift) interval, Ar « 850A - 900A. 
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Lastly, we also expect the mean free path to increase with decreasing redshift as the 
universe expands. In the SDSS analysis at 2; ~ 4, the mean free path was so short that we 
had no sensitivity to redshift evolution in kll from a single stacked spectrum (PWO09). By 
examining stacked spectra over a (narrow) range of redshifts, however, we revealed evidence 
for an increasing A^^, with decreasing redshift. In the WFC3/UVIS-G280 stacked spectrum, 
the flux covered includes opacity from absorbers spanning from 2:912 ~ 1 — 2.5 and one 
must consider explicit redshift evolution in Kll. We assume the iigi2{z) opacity term in 
Equation \T2\ evolves as a power- law with the scale-factor a = 1/(1 + z), 

i^912{z) = '^912(^stack) 

where ^^stack is the median redshift of the quasars whose spectra are averaged. This assumed 
functional form for kgi2 is motivated by the fact that the evolution should be dominated by 
expansion of the Universe and the empirical observation that many aspects of the IGM are 
well-modeled by power-laws of the scale factor. We would reconsider this parameterization 
if it provided a poor model for the data. 



l + z 

1 ~l~ 2^stack 



(17) 



5.2. Generating the WFC3/UVIS-G280 Stacked Spectrum 

As described in the previous subsection, constraints on the mean free path at z < 2.5 
may be derived from analysis of the average quasar spectrum of Zcm ~ 2.5 quasars. Our 
WFC3/UVIS-G280 dataset comprises 53 quasars with Zem = 2.3 — 2.6 (median of 2.44), each 
chosen to be free of significant associated absorption. 

We have taken the following procedure to combine the individual spectra, accounting for 
the differences in Zem and also the non- uniform dispersion of the WFC3/UVIS-G280 grism 
spectra. 

First, we correct for reddening of the quasar flux by adoptin g the E(B—V) estima t e from 



Schlegel et al.l (119981 ) and a standard Galactic extinction curve iFitzpatrick fc Massal (jl990l ). 
In general, this is a small correction because the median reddening value is E{B — V^) = 
0.02 mag. On the other hand, our spectra cover the so-called 2175A bump in the extinction 
curve which implies an additional correction. Second, we generate a rest-frame wavelength 
array with fixed dispersion AA. The dispersion value was set to be large enough to include at 
least one entire pixel from the WFC3/UVIS-G280 spectra at rest wavelengths Ar < 1215A. 
Specifically, we adopted AA = 6.1897A. Third, each quasar spectrum was shifted to the rest- 
frame and normalized by the median flux from A^ = 1450A — 1470A. Fourth, we assigned 
the entire flux of each pixel from the original spectrum to the nearest pixel of the stacked 
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spectrum. Lastly, all of the flux values in the final spectrum were averaged to produce the 
fx spectrum, normalized to unity at A = 1450A. 

The resultant WFC3/UVIS-G280 stacked spectrum is shown as the solid black line in 
Figure [HI and tabulated in Table |6l To assess sample variance, we generated a set of stacked 
spectra with standard bootstrap techniques. Specifically, we generated 500 stacked spectra 
by drawing randomly from the 53 quasar spectra, allowing for duplications. The shaded 
region in Figure |H] shows the RMS at each pixel after clipping any 3a outliers. To some 
degree, this bootstrap analysis provides estimates on the uncertainty in our measurement of 
the average observed quasar flux at z ~ 2.4, subject to the selection criteria of the sample. 

Overplotted o n the stacked spectrum is the average radio-quiet quasar spectrum of 



Telfer et al.l ( 120021 ). produced from a set of 2; ~ 1 quasars observed with UV spectrometers 
on HST and modifled as follows. The Telfer spectrum has been normalized at Ar = 1450A, 
smoothed to the WFC3/UVIS-G280 spectral resolution, and rebinned to the dispersion 
solution of the WFC3/UVIS-G280 stacked spectrum. At wavelengths A^ ^ 1300A, the 
WFC3/UVIS-G280 stacked spectrum has flux that modestly exceeds the Telfer spectrum. 
This suggests that the WFC3/UVIS-G280 stack has a bluer underlying power-law continuum 
but it could also result from strong line-emission (e.g. O I 1302) at those wavelengths. We 
explore this issue in greater depth in § 15.4.21 

At all wavelengths A^ < 1200A, the Telfer spectrum exceeds the WFC3/UVIS-G280 
stack. Here, the WFC3/UVIS-G280 stacked spectrum includes absorption by the intergalac- 
tic medium, whereas the Telfer spectrum has been corrected for IGM absorption (to the 
best of their abilities). For wavelengths Ar ~ llOOA, the offset is ~ 10% and it increases 
with decreasing wavelength as additional terms in the Lyman series contribute (Equation |9]). 
Beyond the Lyman limit, the flux for the WFC3/UVIS-G280 stack drops rapidly owing to 
the integrated continuum opacity of the Lyman limit t^q (Equation [TT]). At Ar ~ 6OOA, the 
ratio of the Telfer spectrum to the WFC3/UVIS-G280 stack implies a total effective optical 
depth tJ^^ ~ 1.5. We expect the opacity at these wavelengths to be dominated by the 
Lyman limit opacity. We now model this absorption and thereby place new constraints on 
the IGM at z ~ 2. 



5.3. Mock Spectra 

To better gauge the effects of sample variance, to test our stacking, and to examine the 
evolution of quantities such as Tg^™^"^, we created a number of sets of mock WFC3/UVIS- 
G280 spectra. The mocks were utilized at each stage of the analysis of the real WFC3/UVIS- 
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G280 stack described below, both to develop and test our analysis algorithms, and to better 
explore the range of likely errors for each measured quantity. 

The mock spectra were g enerated froni simu l ated H I line distributions b ased on em- 
pirical parameterizations (e.g. iMadau et al.lll996t IWorseck fc Pro chaskal 120111 ). Under the 
implicit assumption that the Lya forest is a collection of Voigt profiles with uncorrelated 
parameters (redshift column density A^hi, Doppler parameter b), we populated each sim- 
ulated sightline with absorbers until the Lya effective optical depth converged to a given 
value. If the Lya effective optical depth evolves as rj^" 



oc 



/ oc (1 + ^)^. At ;z < 1.5 we assumed t^^" = 0.017(1 + 2 



1 + z 

1.20 



^"""^ the line density is 
(Kirkman et al. 2007), 
whereas at z > 1.5 we incorporated the observed steepening in the riS'" evolution by taking 



r^^/" = 0.0062 (1 + zf^^ (Dall'Aglio et al. 2008). 

For the Doppler p arameters we adopted the single parameter distrib ution funct i on by 
Hui fc RutleJ fll999l ) dn/db oc b-^exp {-b^ /b^) with b^ = 24 kms"^ kim et all [2OOII ) 
restricted to 10 kms~^ < & < 100 kms~^. The main distribution of interest is the column 
density distribution for which we assumed a triple power law /(A^'hi, A^hi", N^f-^) = CiX A^j^f ' 
over the range 12< logA^Hi <22. The constants Cj implicitly depend on the evolution, 
so we determined them from 5000 simulated sightlines for each of our sets of slopes By 
varying these three slopes we recovered different values of the mean free path (Equation 11). 
In total, we ran 42 different models that resulted in a wide range of MFPs at z = 2.4 from 
40 Mpc to the horizon length. For simplicity we discarded most of these and considered only 
those with simulated MFPs close to our measured value. 

With these Monte-Carlo simulated line lists we then generated mock WFC3/UVIS-G280 
spectra as follows: For e ach of the 53 quasa rs in our sample we first generated a mock quasar 
SED by multiplying the lTelfer et al.l ( 12002| ) composite with a power law fx = (Arest/2500A)" 
with a drawn from a Gaussian distribution with a = 0.2. The mock SED was normalized 
between 1430 and 1470 A. We then randomly drew one of the 5000 sightline realizations, 
simulated the resolved Lyman series and continuum spectrum at < z < Zem, multiplied 
it onto the interpolated SED, convolved it with the WFC3/UVIS-G280 grism line spread 
function and rebinned it to the WFC3/UVIS-G280 grism dispersion solution. Lastly, we 
added Gaussian noise to the spectra to match the quality of our WFC3/UVIS-G280 spectra 
(S/N~ 30 at I350A), taking into account both the shape of the mock quasar spectrum and 
the WFC3/UVIS-G280 sensitivity function. 
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5.4. Modeling the Stacked Spectrum 

In this subsection, we describe our approach to modeling the stacked quasar spectrum 
from the WFC3/UV1S-G280 sample (Figure [8], Table [6]). Our scientific emphasis is to place 
new constraints on the magnitude and evolution of the mean free path A^^p at 2; ~ 2. 
Nevertheless, the data redward of the Lyman limit offer additional constrains on the intrinsic 
SED of 2; ~ 2.5 quasars and also the Lyman series opacity at 2 ~ 2. We consider each of 
these in turn. 



5.4.1. QSO SED 



Central to our experiment is the fact that we use quasar light to probe the foreground 
IGM. To some extent, it is an observational necessity for revealing the nature of this medium. 
Galaxies are too faint for current facilities, especially at Xj- < 2000A, and gamma-ray bursts 
or supernovae are too faint and/or fade too rapidly to generate sufficiently large samples, 
especially with UV spectroscopy. 

As is apparent from previous work on the rest-frame SEP of quasars at UV wavelengths 



(e.g. iTelfer et al.ll2002l : IVanden Berk et al.ll200ll : IScott et al.ll2004l ). these sources do not have 
smoothly varying, intrinsic continua. Although the underlying SED is roughly a power-law 
{fx oc A° with a ^ —1.3 for Ar > 1200A), there are significant and broad emission lines from 
the H I Lyman series and metal-line transitions. Furthermore, individual quasars exhibit a 
range of power-law slopes and varying strength of line-emission. Despite this diversity, it 
is re markable that quasars as a population have a very similar average SED at all redshifts 



(e.g. 



Vanden Berk et al.ll200ll ). Outside the Lya forest, for example, the SED in the UV 
is nearly identical between the z ~ 1 Telfer spectrum and the average spectrum of z ~ 3 



quasars (IVariden Berk et al.ll2001l) . There is only a hint of a harder SED for higher z quasars 
( Telfer et al. 2002h and also m odest differences in the equivalent width s of high-ion i zation 
emission lines flBaldwinl 119771 1 At wavelengths A^ < 1200A, however, iTelfer et al.l (120021 ) 
report a bre ak in the quasar SED power-law to fx oc A^*^'^ that is not evident in other 
analyses (e.g. lScott et al.ll2004j ). This difference could be r elated to redshift e volution and/or 
the average luminosity of the quasars sampled. Recently, IShull et al.l (120121 ) have published 
a composite spectrum from high-dispersion HST/COS observations (with corrected Lyman 
limit absorption) and find results that are in good agreement with the Telfer et al. SED. 

Because of the remarkable similarity in the average SED of bright quasars with redshift, 
we adopt the following assumption for modeling the WFC3/UVIS-G280 stacked spectrum: 
we assume the underlying SED is the Telfer et al. radio-quiet spectrum modulated by a tilted 
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power-law, 



SED 
A 



Telfer 
A 



A 



5aT 



1450A 

and a scaling parameterized by Ct- Returning to Figure |8|, we propose that the offset 
between the Telfer spectrum and the WFC3/UVIS-G280 stack at Ar ~ 1300A may result 
from a difference in the average power-law of the two stacked spectra. Indeed, the data are 
well matched if we assume 6aT ~ —0.5. Interestingly, such a tilt is akin to arguing that 
there is no break in the average quasar SED at A,. ~ 1 200A, contrary to the findings of 
Telfer et al.l ( 12002| ) but consistent with other estimations ( IScott et al.l 12004 ). 



The inclusion of a tilt also allows for the fact that iTelfer et al.l ( 12002| ) may not have 
properly or entirely corrected for IGM absorption at A < 900A. For example, the authors 
did not correct for partial LLS, which are expected to be rare at 2; < 1 ( iRibaudo et al.ll201l[ l 
but not negligible. 

In the following, we restrict 6aT by demanding that the effective Lya opacity t^^" de- 
rived from our st acked spectra matche s previous estimates from analysis of higher spectral- 
resolution data (IKirkman et al.l 120051 ). Similarly, the parameter Ct allows for a normal- 
ization offset between the Telfer stack and the WFC3/UVIS-G280 stack. Although each 
spectrum was normalized at Aj. ~ 1450A, there may be a modest difference in the emission- 
line strength at these wavelengths and it is non-trivial to precisely measure the flux of each 
stacked spectrum at a given wavelength. Therefore, we allow for a 5% modulation in the 
normalization. 



5.4-2. Constraints on the Intrinsic SED from 

As noted in the previous sub-section, we have significant reason to believe that the 
average intrinsic quasar SED for our WFC3/UVIS-G280 sample does not follow the standard 
Telfer quasar spectrum. This is suggested by the offset at Ar ~ 1300A, but is even more 
apparent in the Lya forest (Figure [8]). At Ar ~ llOOA, for example, the Telfer spectrum 
lies only a few percent above the stacked spectrum which would imply almost zero effective 
Lya opacity at z ~ 2. This contradicts our knowledge of the IG M at these redshifts 



as informed by high S/N, high spectral- resolution observations (e.g. iKirkman et al.l 12005 



Kim et al.l 120071 ). Indeed, we now invert the problem to constrain the average quasar SED 



of our quasar sample, using previous estimates of t^". 



Previous authors have estimated the effective opacity from the Lya forest using higher 
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spectral resolution data of z ~ 2.5 quasars (IKirkman et al.l l2005l : iKim et al.l 120071 ). The 
standard approach is to estimate the intrinsic quasar continuum f^^^ from the data directly 
and then measure the average absorption, 



obs / fSED 



1 - rvf 



(19) 



between the Lya and Ly/3 emission lines corresponding to a redshift slightly less then z^m- 
Depending on the author, estimations for Da may include only the low density Lya forest, 
all Lya lines (i.e. including strong absorbers like damped Lya systems), and/or metal- line 
absorption. Because our WFC3/UVIS-G280 stacked spectrum in cludes the opacity f rom all 
of these sources, we wish to compare against a total estimate D^. iTytler et al.l (120041 ) report 

2.14 of 0.1 8 ± 0.04. Analysis of SDSS quasar spectra provides a 



a total value at z 



similar value ( iDairAgho et al.ll2008l ) 



In Figure [9] we present a series of SED models for our intrinsic spectrum comparing 
to a range of tilts and normalizations {Sa^, Ct) applied to the Telfer spectrum. In each 
case, we have demanded that the average opacity at Ar = [1080, 1140]A fall within the la 
interval given by D^. Restricting Ct to ±5%, we find (5aT values ranging from f» —0.2 to 
— 1. It is very unlikely that the unaltered Telfer spectrum provides a good description of 
the intrinsic SED for our WFC3/UVIS-G2803 cohort. Instead, the data favor a bluer SED 
at Ar < 1200A. In fact, the preferred tilt for Ct = 1 of da^ ~ —0.5 nearly co rresponds to 



maintaining the f\ oc A pow er-law that is observed for quasars at X, > 1200A flTelfer et al.l 



2OO2I : IVanden Berk et al.ll200ll ). While it is possible that a portion of this tilt relates to error 
in the fluxing of the WFC3/UVIS-G2803 spectra (Paper I), it is very unlikely to be entirely 
explained by such systematic effects. Instead, we conclude that the WFC3/UVIS-G280 
quasar cohort has a harder SED than the 2; ~ 1, Telfer et al. radio-quiet sample. In addition 
to the implications for our analysis of the IGM, a harder SED would imply a higher emissivity 
of ionizing photons from z ~ 2 quasars. In turn, it would increase the quasar contribution 
to the extragalactic UV background at these redshifts. 

The bluer SED relative to the Telfer analysis may also arise from the color- select ion of 
SDS S quasar candidates. At 2; ~ 2.5, the average quasar color lies near the stellar locus 



[e.g. [Richards et al.ll2006l ) and therefore, the SDSS team weighted their targeting algor ithms 
toward quasars with UV-excess (i.e. bluer SED). Indeed, IWorseck &: Prochaskal (120 111 ) have 
estimated that SDSS quasars with spectroscopic redshift z ^ 2.5 would have Sa^ ~ —0.3 
(see their Figure 16). We encourage additional analysis of quasars at 2; > 2 to further explore 
this issue. 



Although the above analysis prefers a hard SED, we emphasize that a tilt 6aT ~ — 1 is 
not well supported by the individual spectra. None of the quasars exhibit such a blue SED at 
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A < I200A. For example, the quasars J083326+081552 and J121519+424851, which appear 
to have minimal Lyman limit absorption from the IGM, are better described by Sa^ ~ —0.5 
at such wavelengths. Therefore, in the following, we restrict the analysis to SED models 
with 6aT = [-0.8,-0.2]. 



5.4-3. Lyman Series Opacity 

Although the Lya forest is characterized by a series of narrow [Sv < 50km s~^) and 
stochastically distributed absorption-lines, the WFC3/UVIS-G280 stacked spectrum exhibits 
no discrete absorption features. This follows from the low spectral-resolution of the individ- 
ual WFC3/UVIS-G280 spectra and also the effects of spectral stacking. As detailed in § I5.H 
the Lyman series opacity should increase with decreasing rest- wavelength as additional tran- 
sitions contribute (Equation [9]). Below the Lyman limit, however, the total Lyman opacity 
r^™'^'^ is then expected to decline with decreasing redshift (i.e. as the universe expands). 

One may relate the effective opacity at a given Lyman transition to the frequency 
distribution of IGM absorption lines f{Nm, b, z) as. 




^eff= / / f{Nui,b,z)exp{-T:)dNnidbdz , (20) 



where r" is the line opacity of transition n and is a function of iVni, b, and z. Because 
f{Nm) is estimated to decline at least as steeply as N^^'^, one predicts the effective opac- 
ity of a given Lyman transition to be dominated by the the minimum column density for 
line saturation. This means that rjj^" is set by lines with Nm ^ 10^'^ cm~^ whereas the 
higher order transitions are defined by lines with A^hi ~ lO^^cm"^. Indeed, for reasonable 
f{Ni{i,X) distributions we estimate that roughly half of the total opacity derives from the 
first few transitions (tt. = 2 — 5, i.e. Lya-Lj6) with the remainder contributed by higher order 
transitions. 

To estimate the t-^™*^"^ opacity from the frequency distribution, therefore, one requires: 
an accurate description over many orders of magnitude in Ahi, an assessment of the Doppler 
parameter distribution, and also an accurate estimation of their evolution with redshift. 
While there has been some analysis on f{N, b, z) at z ~ 2, the constraints are quite limited. 
Given this uncertainty, we have made our own estimate of t^™^" from the WFC3/UVIS- 
G280 stacked spectrum, corresponding to 2; = ^stack = 2.44 . In Figure [TO] we estimate t^™*^"^ 
allowing for the range of intrinsic quasar SEDs described in the previous subsection. The top 
panel shows the recovered t^™"^" values as a function of tilt Sa^ in the SED. The bottom 
panel shows the range of allowed SEDs compared against the WFC3/UVIS-G280 stacked 
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spectrum. We also present the unaltered Telfer spectrum (purple, dotted line) to further 
emphasize that this SED is inconsistent with the observations. If we assume that Sa^ = —0.8 
to —0.2 at 95% C.I., then this implies t^™^° = 0.40 ± 0.15 at similar confidence. There will 
be additional uncertainty from sample variance, but we estimate that this contributes less 
than the uncertainty related to the quasar SED. 

This estimate for t^^™^" may be crudely compared against the incidence of r > 1 
LLS derived from our spectra. Specifically, one may assume that ~ 50% of the opacity is 
contributed by the higher order lines and then compare this value to an estimate for r^™'^'^ 
for LLS from an evaluation of equation [TTl Taking A^hi > 10^'''^^ cm~^, i{z) = 1.5 for z = 2.0, 
bes = 35 km s~\ and assuming f{N, z) oc N~^-^{l + zY-^'^ yields Teg = 0.06 for the transitions 
= 5 — 30. This value is considerably lower than the t^^™^" value derived from above, and 
it suggests that systems with A^hi < 10^^'^ cm^^ dominate the Lyman series opacity. This 
conclusion is c onsistent with our in ferences of a steep f{N) distribution at A^hi 
see also lRibaudo et al.ll201ll ). 



10^^ cm-2 



In the following, we model the WFC3/UVIS-G280 stacked spectrum with a range of 
T^l^^^ values. In fact, we take exactly the r^™^" value required to reproduce the observed 
flux at Ar = 912A for a given quasar SED and then allow for a 10% variation which signifi- 
cantly exceeds what can be measured from our high S/N stack and therefore is a conservative 
range. For < 912A, we further assume that t^J^™"^" decreases as a (1 + z)'^^ power-law. 
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(21) 



This is justified by the fact that /{Nm, z) is generally well-described by such a power-law. 
On the other hand, di fferent A^ht regimes a nd different redshift ranges may be described 
by different 7,- values (jjanknecht et al. 2006 ). We proceed under the expectation that the 
total T^™^" is nevertheless well described by a power-law over the redshift interval relevant 
to our analysis 2912 ~ 1.5 to 2.4. Specifically, we allo w 7t to range from 1. 2 to 1.8 which 
spans the range of v alues estimated for the Lya forest (jJanknecht et al.ll2006l ) and LLS (§12 



Ribaudo et al. 



201 If ) at 2 < 2. 



5.4-4- Lyman Limit Opacity 

Following the formalism presented in § 15. Ij we characterize the Lyman limit opacity kll 
as given by equations [T2] and [T71 This is a two parameter model that sets the normalization at 
z = -Zstack and allows for redshift evolution. In the following, we will demand that 7^ > 0, i.e. 
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that the A^\^ increasecl with decreasing redshift as expected from the decreasing incidence 
of LLS (§11X1]). 

5.5. Constraining the Mean Free Path 

We now proceed to compare a suite of models against the observed WFC3/UVIS-G280 
stacked spectrum to estimate the mean free path A^jp. We proceed with standard analysis 
where the uncertainty at each pixel in the stacked spectrum is estimated from a bootstrap 
analysis (§ I5.2p . This yields a 'best' estimate for A^^p. In turn, we provide an estimate 
of the uncertainty in this quantity by repeating such analysis on a set of stacked spectra 
generated with standard bootstrap techniques. There are several issues to note regarding 
this approach. First, the bootstrap analysis of the stacked spectrum provides an estimate 
of the RMS at each pixel in the stack, but we caution that the PDF need not follow a 
true Gaussian. On the other hand, we find our results are relatively insensitive to what 
we assume for the error in the stacked spectrum. Second, the bootstrap analysis may not 
properly reflect the full uncertainty in the stacked spectrum related to sample variance, i.e., 
the stacked spectrum was derived from a total of 53 quasars. Lastly, this approach ignores the 
fact that the data at Ar < 9 12 A are highly correlated. This correlation occurs because of the 
nature of Lyman limit absorption; a system a.t z = Zgi2 attenuates the flux at all wavelengths 
Ar < 912A(1 + Zgi2)/{1 + Zem)- Therefore, one predicts (and observes) a monotonic decrease 
in the stacked spectrum for Ar < 9 12 A. The first two issues described above suggest an 
underestimate of the uncertainty while the last point may lead to an overestimate. 

Our model has six parameters: two for the quasar SED {Sa^, Ct), two to model the 
Lyman series opacity {tIq'^^'^{z = -Zem), 7t), and two to model the Lyman limit opacity 
('^9i2(^stack), 7k). In the previous subsection we imposed constraints on the parameters as 
summarized in Table [3 We then constructed a g'^id in this six dimensional parameter 
space for rest wavelengths Ar = 700.0 — 911. 76A. Although the stacked spectrum extends 
to Ar = 600A, uncertainty in the quasar SED, the evolution of t^'^^^^, and sample variance 
are much greater at these wavelengths. Furthermore, our principle goal is to estimate A^^, 
which we find occurs at Ar ~ 800A. These issues motivated our decision to terminate the 
quantitative comparison at Ar = 700A. 

Figure [TT] shows the results for the best-fit model and the suite of models with Ax^ = 
xl — xlmin < This choice for a limiting Ax^ is primarily illustrative; it does not impact 
our uncertainty estimate for A^£p. The top panel displays the range of allowed 7^ and 



Formally, 7^ = also implies an increasing because of the expanding universe. 
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'«9i2(-2stack) values; there is an obvious degeneracy between these two parameters. In short, 
the data permit a more opaque universe ai z = 2.4 that rapidly evolves to a lower opacity 
or a less opaque universe that evolves more slowly. 

Despite this degeneracy in «:9i2(-2stack) and 7^, the results imply a relatively tight con- 
straint on A^fp. In Figure [121 we plot the Zgf^^ values for the full set of allowed models 
as a function of the SED tilt Sa^ and redshift evolution of the Lyman series opacity 7,-. 
We find that Zq^2 values from 1.8 - 2.0 which for our adopted ACDM cosmology implies 
■^mfp — ~ 300/172^ Mpc. Again, we consider these estimates to be largely illustrative, not 
quantitative; systematic effects and sample variance dominate the uncertainty. 

One approach to assessing uncertainty related to sample variance is to repeat the above 
analysis for a suite of stacked spectra generated by standard bootstrap techniques. Fig- 
ure [12] presents the A^fp values for each of 500 realizations of the stacked spectrum us- 
ing the same formalism described above. We recover a median (mean) mean free path 
of A^fp = 243 (252) /i72^ Mpc. We consider this value to be our best estimate for A^^, at 
z = 2.44. This same analysis suggests an uncertainty in A^fp of ~ 20%. This should be 
considered a minimum estimate for the uncertainty. We also considered the systematic er- 
ror reddening corrections might have on our stack analysis by turning off the de-reddening 
and re-computing A^£p. The results of this change are within our estimates on the error for 
A^£p, with the implied value approximately 15 Mpc lower than our fiducial value. The set of 
best-fit parameters are presented in Table [3 

Consider, further, the sources of error contributing to the X^^^ measurement. Our 
analysis only considers the stacked WFC3/UVIS-G280 spectrum at rest wavelengths Ar = 
700 — 910A where one observes a rapid (> 60%) decline in the relative flux with decreasing 
Ar. Our model of the relative flux consists of three independent effects: (i) a reduction in 
the flux due to the effective opacity of Lyman limit absorption r^; (ii) an increase in the 
flux due to decreasing Lyman series opacity t^™*^"^; and (iii) changes in the flux owing to 
relative changes in the QSO SED. We emphasize that the latter two effects are expected to 
be small (< 10% in relative flux), in particular because the analysis is performed on such 
a narrow range of A^.. For example, if we ignore variations in t^™^", we recover the same 
A^£p value to within 3%. Similarly, varying the tilt in the QSO SED modifies A^^^ only 
several percent. These two aspects of the model, however, do work together in that they 
both give higher relative fluxes at lower Ar. We find that adopting a redder QSO SED (e.g. 
the Telfer spectrum without any tilt) and no 7-^™'^'^ evolution would yield an f» 10% higher 
A^fp value than our favored value. Even redder (i.e. 6aT> 0) tilts have been considered, and 
would significantly increase our estimate of A^fp. These tilts are highly disfavored, however, 
as they summarily fail to reproduce the constraints provided by r^°. We cannot reject the 
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possibility that the underlying QSO SED could have additional strong inflections near rest 
frame wavelengths of 912A, thus adding significant additional uncertainties to our estimate 



of X^ f „, but we see no indications of su ch features in other QSO SED studies (ITelfer et al. 



— 

2002; 


Scott et al. 


2004; 


ShuU et al. 


2012) 



aspects of the modeling impose an f» 10% systematic error in the X^l analysis. 



Presently, uncertainty in the estimated X^^^ value is dominated by sample variation. 
As indicated in Figure [HI the stacked spectrum has an ^ 20% scatter at ~ 800A. This 
dominates the uncertainty in the results (Figure US]). We estimate that one would have to 
increase the sample size by at least a factor of four before uncertainties related to the QSO 
SED or r^™'^" would contribute substantially. Lastly, we comment that extending the A^fp 
analysis to lower redshifts {z < 1.5) where the A^^p value is presumably much larger will 
prove progressively more difficult. As the Lyman limit opacity has a weaker impact on the 
relative flux, uncertainties in the QSO SED and/or r^™*^"^ could dominate the analysis. 



5.6. Comparisons with other estimates of X^f^ at 2; ~ 2.4 



Our bootstrap estimate of A^£p= 252 ± 52hj2 Mpc can now be compared with other 
estimates derived through different methods at this redshift, namely through application 
of Equation [TT] after adopting an estimate for f{Nm,X)The exact approach to incorporat- 
ing systems with ^ 1 can dramatically affec t the results, and lend to ver y different 
estimates of A^fp when using this methodology. iFaucher-Giguere et al.l (l2008al ) arrive at 
a value of A^^fp= 163/;,f2^ Mpc for z = 2.4, but establish the nor r aaliza tion of f(Nui,X) 
by considering only those LLS with > 1. ISongaila fc Cowid (120 lOh provide new ob- 
servations and analyses of LLS at z > 4, along with low redshift GALEX LLS observa- 
tions to adopt a single power-law form for / (A^ht,^) in the LLS regime, and the same 
methodology as iFaucher-Giguere et al.l (l2008al ) to arrive at X^^l= 172/172^ Mpc at ;z = 2.4. 
Worseck fc Prochaskal ( 120111 ) employ the A^fp from PWO09, plus new constraints on /(Ahi, X) 
at 14.5 <logA'Hi< 19.0 to determine a multiple power -law f(Ny['[,X), and arr ive at a 
higher value of A^^^ ^ 220hj2 Mpc at z = 2.4. Finally, Iflaardt fc Madaul J2012h employ 
the f{Niii,X) from POWIO to arrive at X^J^l= 185hji Mpc at z = 2.4, with a very steep 
redshift dependence at these redshifts (by z = 2, A^fphas increased to 242h^2 Mpc). Clearly, 
the derived A^fp values from /(Ahi,-^) estimations are very sensitive to the poorly con- 
strained regime at r^^^ ^ 1. The methodology presented here has the advantage that it is 
wholly independent of /(A^hi, X) and the results for A^fp can then be used to offer constraints 
on /(Ahi,-^), as in the next section. 
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6. Constraints on f{N) at 2; ~ 2.2 

We now turn our focus to an exploration of tlie column density distribution function 
f{N}ii,X). Our analysis of i{X) for the LLS along with our stacked spectrum analysis of 
A^fp allow us to place constraints on f{Niii,X) over many decades in H I column density. 
We have already seen in section 14.1.11 that we expect significant (and multiple) deviations 
from a single power-law for f{Niii,X) in order to match constraints from the Lya forest on 
the low H I column density end, and the SLLS and DLA on the high H I column density 
column end. We address the problem in stages, from high to low A^hi before performing a 
full and simultaneous fit to all of the observational constraints at 2; ~ 2. 



6.1. Constraints on f{Nm,X) for Nm > lO^^cm-^ 

Our survey of Lyman limit systems provides an integral constraint on the A^hi fre- 
quency distribution at high values (equation HI POWIO). Specifically, we observe an in- 
cidence of r > 2 LLS at 2 ~ 2 of ^(X)t->2 = 0.29 ± 0.05 corresponding to systems 
with Ahi > lO^^'^cm"^. This integral constraint may then be compared against the ob- 
served incidence of super Lyman limit systems (SLLS; 10^^ cm~^ < A^hi < 10^°'^ cm"^) 
and damped Lya systems (DLAs; A^hi > lO^"''^ cm~^) that have been est i mated from sur- 



veys analyzing the dampirig wing of the Lya transitio n (jPeroux et al.ll2005t iProchaska et al. 



2005; O'Meara et al. 2007 



Prochaska fc Wolfd l2009l ) . If one assumes a functional form for 



f{Niii,X) at Ahi < 10^^ cm ^, e.g. a single power-law, then we may directly constrain its 
parameters as follows. 

Figure [7] presents our estimates for £{X)r>2 compared against the incidence of SLLS 
and DLAs, in c umulative form. For the D LAs we have assumed ^(Ar)DLA = 0.05 from the 
SDSS survey of IProchaska fc Wolfd (120091 ). derived from the redshift interval z = [2.2,2.4]. 
We note that s ystematic error associate d with the SDSS spectral coverage affects iuLAi^) at 
these redshifts (iNoterdaeme et al.ll2009l ). but the DLAs have a sufficiently small contribution 
that we may neglect this issue. We also adopt the /(Ahi,^) distribution for DLAs from 
Prochaska &: Wolfd (120091 ) over the same redshift interval. This is plotted in Fi gure [TH For 
the SLLS, we adopt the value of ^(X)slls = 0.13 given by (lO'Meara et al.ll2007l ). which they 
derive for SLLS with log A'hi> 19.0 and 1.7 < ;z < 3. 

Together, Figure [7] shows that the SLLS and DLAs contribute 50 — 80% of i{X)r>2- 
Even if we assume the upper end of our £(X)^>2 estimate!^ it is evident that systems with 



10-1 



We further note that lRibaudo et al.1 (|201lh reported e{X)r>2 = 0.27±0.09 from their survey of the HST 
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A^Hi = lO^^'^cm"^ — lO^^cm"^ have a relatively modest contribution. In turn, this implies a 
shallow /(A^Hi; X) for ^hi < 10^^ cm~^. Figure shows the range of constraints for a single 
power-law covering that A^hi interval, 

f{10^'-' cm-2 < Nni < lO'' cm'^; X) = ^lls^^"' , (22) 

which assumes (i) £{X)r>2 = 0.29±0.05; (ii) W(X)=0.05; (iii) log/(ArHi = lO^^cm-^X) = 
-20.2±0.2. The latter constraint follows from £(X)slls = 0.13±0.04 and /3slls = -1.2±0.2. 
We estimate log/cLLs = —9.2 and /3lls = ~0.9;1^q;q5, and note that the two parameters are 
highly correlated. 

The derived slope for f{Niii,X) in the LLS regime is very shallow; /3lls > ~1 implies 
the universe exhibits greater cross-section to gas with A^hi = 10^^ cm~^ than 10^^ cm~^, which 
seems unlikely. However, a similar result was derived (in nearly identical fashion) for LLS 
at z ~ 4 by POWIO. At both epochs, the relatively low incidence of r > 2 LLS implies a 
flattening of /(A^hi, X) at A^hi ~ 10^® cm"^ that we speculate is associated with transitioning 
from an optically thick to an optically thin regime. We return to this issue in § [71 

6.2. Constraints on f{NHhX) for A^^hi ~ 10^^ cm'^ 

In Section H] we presented results on the incidence of LLS for several limiting opacities 
at the Lyman limit: rnmit = 0.5, 1 and 2. While these measurements are not independent^ 
(they were derived from the same sightlines and the samples overlap), one may compare the 
results to offer an estimate for f{Ni{i,X) in the interval A^hi = lO^^'^cm"^ — 10^^'^ cm~^ 
corresponding to = 0.5 — 2. 

Once again, we parameterize f{Ni{i,X) as a single power-law, 

filO''-' cm-2 < AThi < lO^^-^ cm-^ X) = fcpLLs<r' • (23) 

To estimate both /3pLLS and /cpLLS, we must impose an additional constraint on /(A''hi,X) 
from the LLS analysis of the previous sub-section. Specifically, we consider the range of 
allowed values for /(Xhi = 10^^'^ cm~^, X) which Figure [H] demonstrates could range from 
^ 10~^^'^ — 10~^° cm^. In Figure [TBk. we plot the offset in incidence Ai{X) from the th^h > 2 



archive. 

^^We further remind the reader that £{X)t->o.5 was required to exceed £{X)r>2 in the analysis. While this 
must be true, it further emphasizes that the two results are highly correlated. 
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measurement for rnmit > 1 and rnmit > 0.5. Overplotted on these offsets are a series of curves 
for the predicted offsets, as a function of rumit, for a range of /(iVni = 10^^'^ cm~^, X) and 
/3pLLS values. Models with /(A^m = 10^'^'^ cm"^, X) < 10"^^ cm^ predict too few pLLS for 
/5pLLS > ~2. Indeed, only for /(Xm = 10^^-^ cm^^, X) = 10^^® cm^ is a power-law with 
/3pLLS ~ — 1 permitted in the pLLS regime. 

This issue is further illustrated in Figure [T5b where we indicate the allowed values 
for /SpLLS (i-e. in agreement within la of both A£{X) values) as a function of /(Xhi = 
10^^'^ cm~^, X). We conclude that /plls(^hi, must be steeper than /3pLLS = —1-5 unless 
f{Nm = 10^^-^ cm"^, X) > 1 0"^^ '' cm^. A s t eeperi ing of /(Xhi,X) in the pLLS regime has 



been previously reported by iRibaudo et al.l (1201 ll ) based on their LLS analysis. Similarly, 
POWIO argued for a steepening of /(Xhi,X) at Xhi < 10^''cm~^ in order reproduce the 
mean free path estimates of PWO09. We discuss this result further in § [71 

While a steep /(Xhi,X) is permitted physically, it is quite surprising given the very 
shallow /(Xhi,X) estimated for the LLS systems in the previous sub-section. Indeed, at its 
extreme, the results imply a shift in slope A/3lls = /^lls ~ /^pLLS of unity or greater over 
an Xhi interval of ~ 1 dex. This point is further illustrated in Figure [12] where we estimate 
A/3lls taking into account the i{X) measurements and constraints on the SLLS (the DLA 
results are inconsequential). For our best estimate of £(X)t->2=0.3, we recover A/3lls > 1, 
i.e. a rapid steepening in /(Xhi, X). Only at the upper end of permissible values for i{X)r>2 
may A/9lls be small. We further emphasize that any new measurements of £slls(-^) or 
^dla(-^) that gave higher values would further increase A/^lls- 



6.3. Constraints on /(Xhi,X) for N^i < lO^^cm-^ 

We complete our piece-meal examination of /(Xhi, X) by imposing constraints from the 
mean free path analysis of § [5] Our evaluation of A^£p yields an additional integral constraint 
on /(Xhi,X) as described in equation[TTJ Unlike i{X), the A^jp constraint includes redshift 
evolution in the Xhi frequency distribution. Our results on ijX) suggest very weak redshift 



evolution for /(Xhi,X) at z < 2.5 (Figure [6l see also iRibaudo et al.l [20 111 ), at least for 
systems with Xhi ~ 10^^ cm^^. Furthermore, our estimate of the mean free path is sufficiently 
short that we are only sensitive to /(Xhi,X) for a small redshift interval aX z ^ 2. We 
proceed, therefore, by assuming no evolution in /(Xhi,X). 

Figure [17] shows calculations of the effective Lyman limit opacity for a series of 
/(Xhi,X) models, identical for Xhi > 10^'''^cm~^ but with strict power-laws at lower Xhi 
with /SpLLS = [—1, —1.5, —2, —2.5]. In each case, we estimate from z = 1.86 to z = 2.44 
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which is the redshift interval for which our mean free path analysis yielded t^^ = 1. The 
red shaded region indicates the la uncertainty in this estimate. 

For the favored values of A^jp and for the f{Niii,X) estimates at A^hi > 10^'''^cm~^, 
we find that LLSs with r > 2 contribute approximately half of the required optical depth. 
Indeed, unless we adopt the largest A^fp allowed by the data and the highest incidence of LLS 
then systems with A^hi < 10^ ''•^ cm~^ must make a significant contribution to r^. Figure ITTl 
indicates that /3pLLs ^ in fact /SpLLS < —2 is preferred for our central values. This offers 
(nearly) independent evidence that f{Niii,X) steepens significantly in the pLLS regime. 



6.4. Joint Constraints on a Complete f{Niii,X) Model 

As the previous subsections reveal, our estimates of the mean free path A^fp and the 
incidence of LLS i{X) provide new constraints on the A^hi frequency distribution f{Niii,X) 
at 2; ~ 2.5, especially for A^hi ~ 10^^ — lO^^cm"^. By combining these measurements 
with other evaluations and constraints on f{Ni{i,X) from the literature, we may infer the 
A^Hi distribution across ~ 10 orders of magnitude. This eval uation informs astrophy sical 



quantities sensitive to properties of the IGM (e.g. the EUVB; iHaardt fc Madaull2012l ) and 
provides insight into the physical origin of absorption line systems. At the most basic level, 
we test whether a single f{Ni{i,X) model can reproduce all of the observational constraints. 
Any significant inconsistency would stress a fundamental flaw in our standard methods for 
studying the IGM. 

Before proceeding, it is important to emphasize what may be a subtle distinction be- 
tween two conceptualizations of f{Niii,X): (1) an analysis of the f{Niii,X) distribution 
that one uses to estimate quantities like rj^° and A^£p (Equations [9l[TT|) and (2) an obser- 
vational f{Ni{i,X) derived through line-fitting analysis of absorption-line spectra. At larger 
A^Hi values, where the systems are rare and line-blending is uncommon, the two types of 
f{Niii, X) should equate. At low A^hi, however, the two may diverge. Consider, for example, 
a mock spectrum generated from an analytic /(A^hi? X) following the procedures described in 
Section 15. 3[ One could then perform a line-profile analysis of the mock spectrum and would 
certainty find fewer low A^hi lines than inputted. Ultimately, neither type of /(A^hi, X) actu- 
ally follows the physical description of the IGM suggested by cosmological simulations, i.e. 
the concept of an undulating density field. In this respect, any estimate of f{N}ii,X) has 
limited physical significance. We proceed, nevertheless, with the primary goal of estimating 
an analytic-oriented f{Ni{i,X) distribution. 



Absent a physically motivated model for /(A^hi, X), we assume that it follows a series of 



- 31 - 



power-laws, monotonically decreasing in value with increasing A^hi but discontinuous in the 
first derivative. This approach is partly motivated by previous estimations of /(A^hi? ^) over 
modest A^hi intervals, but it primarily reflects our ignorance of the underlying distribution 
function. Analytically, a power-law approach has the advantage of being robust to non- 
physical deviations, e.g., significant wiggles that a spline or higher-order polynomial may 
introduce. We normalize f{N}ii,X) at A'hi = 10^^ cm~^ with the value lO^^^^. We do not 
evaluate f{Niii,X) for A^hi < 10^^ cm~^ because there are very poor constraints at these 
column densities; for example, gas with such low column densities should contribute less 
than 5% to the total effective Lya opacity!^ We then describe f{Niii,X) as a series of 
broken power- laws with slope /3 in each segment defined by a series of A^hi 'pivots'. For 
example, P20.3 defines the slope from A^hi = 10^°'^cm~^ to the next A^hi pivot. Initially, 
we also allowed for redshift evolution in f{Niii,X), in particular to model A^fp. We found, 
however, that the data offered only very weak constraints. In the following, we assume a 
(1 -|- z)^'^ evolution in f (N}^-\,X), which is ro ughly consistent with the observed evolution in 



at these redshifts (jKirkman et al.ll2005l ). 



In addition to the observational constraints imposed by the results of this paper (A^£p, 
£(X)), we also in clude (1) evaluations of f{Niii,X) estimated from the line-fitting survey of 
Kim et al.l (120021 ) but re-evaluated to our own choice of A'hi binning and to our assumed cos- 



mology. We restrict the evaluations to A^hi < 10^^'^ c m ^ to maintain significant statistical 



power; (2) measurements of f(Nuu X) for SLLS from O'Meara et al.l (120071 ): (3) measure- 
ments of f(Niii,X) fo r DLAs from Prochaska fc Wolfe ( 20091 ): and (4) an estimate of r^J^" 
(from Da) at z = 2.4 (jKirkman et al.ll2005l ). Table O lists these constraints and we discuss 
a few modifications to the reported values and errors below. 

From the empirical constraints, we construct deviates between observation y and model 
m, d = {y — m)/a, with a the 1-sigma uncertainty. This includes the integral constraints 
on f{N}ii,X): t^", i{X) and A^\p. For all of these, we assume Gaussian errors which is 



not fo rmally correct but offers a fair approximation. For the case of t^^", iKirkman et al. 



(I2OO5I ) estimated an error of only 3% but we adopt a more conservative error of 10% to 
account for systematic uncertainty and the fact that our model assumes a fixed Doppler 
para meter for the Lya fores t of 6 = 24km s~^. We adopt a model with the same constraint 
as the lKirkman et al.l (120051 ) analysis that was restricted to systems with A^hi < 10^^'^ cm~^, 
but which includes metal line absorption. For the highest A'hi values (constrained by the 
DLAs), we calculate the deviate from A^hi = 10^^'^ — 10^^ '^ cm~^ and consider this as a single 
constraint. Lastly, we optimized our parametric model with the software package MPFIT 



^^Of course, such gas may trace the majority of the volume of the universe and could contribute a significant 
fraction of the mass density. 
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with initial guesses estimated through a by-eye comparison to the various constraints. 

It is illustrative to consider a series of /(A^hi? X) models with increasing complexity, i.e. 
a increasing number of power-law segments, to establish the simplest model which satisfies 
all constraints. Following the results of the previous sub-section, we model /(A^hi,-^) at 
A^Hi > lO^^cm"^ with a series of 3 power-laws, defined by A^hi pivots at logA^ni = 20.3 
and 21.5 dex. We then optimized a 4-parameter model for f{N}ii,X) under the full set of 
observational constraints. This model assumes a single power-law with slope P12 from A^hi = 
10^^ — lO^'^'^cm"^. Not surprisingly, the model yields very poor results with a reduced chi- 
squared xl > 10. The model cannot match the shallow slope suggested by i{X) (Figure [H]) 
with the steep slope implied by A^^fp (Figured?]) and the Lya forest observations. In short, 
the data require a break in f{N}ii,X) below A^hi = lO^'^cm"^. 

We next evaluated a set of 5-parameter models (4 power-laws) with an additional A^hi 
pivot at ~ lO^^'^cm"^. To our surprise, we found this model well-reproduced {xl ~ 1) 
the Tg'f", X^J^l, i{X), and SLLS/DLA constraints provided /3i2 ^ -1.65 (Figure [Si). This 
implies a significant break to a shallower f{Niii,X) (A/3 > 0.6) for the LLS. Perhaps coinci- 
dentally, this slope is in good agreement with the sh ape of f{Niii,X) reported for the z <1 
universe (jPenton et al. 2004 : Lehner &: Howk 2007 ). This model also implies that absorp- 



tion with low A'hi values (< 10^^'^ cm ^) contribute significantly (~ 50%) to the effective 
Lya opacity. In turn, the model severely over-predic t s the incidence of such weak absorbers 



compared to the empirical estimations of iKim et al.l (120021 ). Those authors reported a more 
shallow f{N}ii,X) at A^hi < lO^^cm"^, inconsistent with our 5-parameter model. On the 
other hand, such low A'hi absorption is the most subject to uncertainties associated with 
S/N, continuum placement, and line-blending. We also posit that this difference may be a 
manifestation of the distinction between the analytic and observational f{Niii,X) distribu- 
tions noted above. We proceed to consider a yet more complex model but note that one is 
not entirely required. 

In an effort to reproduce all of the observational constraints given in Table [9l we exam- 
ined a 6-parameter model (5 power-laws) that includes two A^hi pivots between log A^hi = 14 
and 19 dex. We experimented with parameterizing these two A^hi pivots and found very 
weak constraints on their values. Therefore, we have fixed them at A^hi = 10^^'^ cm~^ and 
A^Hi = 10^ ''•^ cm~^ and caution that the results should not be considered unique. Figure [T8b 
presents the best-fit f{Ni{i,X) distribution compared against the observational constraints. 
The reduced is near unity and is dominated by the /(A^hi,^) evaluations at low and 
high A^Hi- This model shows a shallow distrib ution functi o n at low A^hi which steepens to 



/9i4.5 ~ —2 at A^Hi ^ 10 cm , as inferred by iKim et al.l (120021 ). Again, we conclude that 



f{N}ii,X) is steeper than /3 = —1.6 at A^hi ~ 10 cm and then transitions to /3 ~ — 1 in 
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the LLS regime. This is consistent with our inferences based on the A^£p and partial LLS 
analyses (Sections 15. 5116. 3p . 

Despite the relatively low xt value for this 6-parameter model, there is significant ten- 
sion. A shallow f{N}ii,X) at low A^hi values requires that a significant fraction of the Lya 
opacity arise from systems with A^hi ~ lO^^cm"^. Indeed, a 6-parameter model that pivots 
at Ahi = 10^^ cm~^ instead of lO^'^'^ cm~^ under-predicts the value and n a ay be ruled 
out. The model shown in Figure [TSl meanwhile, over-predicts the iKim et al.l (120021 ) esti- 



mates of f{N}ii,X) at A^Hi > lO^^cm"^. As the precision in the observations improves, we 
speculate that it may become impossible to generate an /(A'hi,-^) distribution that satisfies 
all of the constraints. 

To summarize our main findings on f{N}ii,X) at z ~ 2.4: 
• All 4-parameter (3 power-law) models for f{N}ii,X) are ruled out by the observations. 



A 5-parameter model with a single power-law (/3 



-1.65) spanning from A^hi 



10^2 cm~2 to 



10 



17.5 



cm 



well reproduces the r,g^", A^\^^ £(X), and SLLS/DLA ob- 
servations. This model predicts a large contribution to the effective Lya opacity from 
gas with low A^hi val ues and also over -predicts the incidence of such gas compared 
to the line-analysis of iKim et al.l (120021 ). Such line-analysis is challenged, however, by 
S/N, continuum placement, and line-blending and we recommend further consideration 
of the systematic uncertainties. 



To fit the /(A^Hi; X) evaluations from lKim et al.l (120021 ). one must consider a 6-parameter 

-1.3) at Nni < 



model (5 power-laws). This allows for a shallow f{Ni{i,X) (/3i2 
10^^ cm~^ that must break sharply to a steep distribution (/3i4.5 ~ 
ate Ahi values. 



-2) for intermedi- 



All of the successful models require a relatively steep (/3 < —1.6) distribution function 
for A^Hi = 10^^ — lO^^cm"^ that then flattens in the LLS regime to /3 ~ —1. 



Table [To] lists the best-fit values, the la uncertainties derived from the diagonal elements 
of the covariance matrix, and the correlation matrix for the model. The parameters are well 
constrained, with the exception of the slope at A'ht > lO^^'^cm"^ which is best considered 



as an upper limit (/321.5 < ~3; iProchaska et al.ll2005l ). On the other hand, several of the 



parameters are highly correlated (e.g. ki2 and (3i2) and one cannot assert that our model is 
unique. We also caution that our analysis ignored the parametric freedom of pivot place- 
ment; their values were motivated by the observations. Lastly, we re-emphasize that the 
broken power-law model is assuredly non-physical and therefore an incomplete description 



of the universe. Nevertheless, the results impose a general description of f{Ni{i,X), whose 
implications are discussed below. 



7. Discussion 

We now consider the larger implications of our results. 

7.1. Evolution in i{X) and A|^\^ 

In Figure [191 "we present our measurement of A^jp along with the results of PWO09 at 
higher redshift. Overp lotted on Figure [19] are curves which describe the evolution in A^fp 
with redshift when assume the 5-parameter model functional form for the IGM at z = 2.4. 
The different curves correspond to allowing i{z) to evolve as (1 + 2;)'^, and letting 7 vary 
from 1.5 to 2.5. The curves are forced to have the same value, and match the PWO09 SDSS 
result at z = A. The importance of our new measurement for X^f^ is immediately obvious: 
without information at lower redshifts, the high-redshift data has no discriminatory value in 
7. With the inclusion of the HST data, however, we can rule out extrema in 7, especially 
towards the high end. Specific ally, we can ru l e out 7 = 2.5. We note that our results are 



again consistent with those of iRibaudo et al.l (120111 ). who favor a value of 7 = 1.8. Their 



analyses, however, only explicitly determines the A^^p for systems with r9i2 > 2. Figure 
also shows the need for further exploration of A^jp at other redshifts, specifically at 2; < 2, 
2; ~ 3, and z > 4.5 to provide a full description of the evolution of A^^, over cosmic time. 
Songaila & Cowie (2010) present a measurement of A^fp at 5 < 2; < 6 from ESI spectra 
which indicate a slightly lower value for A^jp than would be predicted from lower redshift 
power-law fits, but we caution that the number of spectra in their sample is relatively low, 
and that their method for determining A^^p involves a number of assumptions not made in 
our analysis here. 

The solid gray curve in Figure [19] shows the Horizon of the universe as a function of 
redshift. When Aj^^p exceeds this radius, ionizing photons are capable of travelling the visible 
universe before being attenuated. In this respect, every source of ionizing photons can 'see ' 



the other. This epoch is termed the 'breakthrough' redshift (e.g. iMeiksin fc Madaijll993[ ). 
Given the current constraints on A^fp and the simple models presented in Figure [191 our best 
estimate is ^break = 1-6 and we require 2;break < 2. Future observations on the incidence of 
optically thick absorption at 2; < 1 should firmly establish the value (Howk et al., in prep.). 
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7.2. The Inflection in f{Nu_i,X) through the LLS Regime 



One of the primary results of our analysis is that the frequency of systems with A^hi ^ 
10^'' cm~^ is significantly higher than that predicted from a simple extrapolation of the 
/(iVHi,X) distribution estimated for A^hi > lO^'^'cm"^. Stated another way, our results 
demand a higher frequency of pLLS than one would predict based on the incidence of LLS. 
We have reached this conclusion from two complementary analysis of our dataset: (i) the 
direct counting of pLLS versus LLS (Figure [TSl): and (ii) our analysis for the total effective 
Lyman limit opacity of LLS (Figure [T71) . Each of these results imply the steepening of 
f{Niii,X) with decreasing iVni. Such a steepening of f{Niii,X) at A^hi ^ 10^^'^ has been 
inferred previou sly at z = 3.7 (P OWIO) based on X^^l^ analysis and at 2; ~ L5 from a direct 
survey of pLLS (IRibaudo et al.ll201ll ). 



This result is further illustrated in Figure [20] which plots the differential contribution of 
systems to i{X) per A log N^i = 0.5 dex interval, Ai{X) /A log A^hi- These Ai{X) /A log A^m 
values were simply calculated from the power-law description of f{N}ii,X) shown in Fig- 
ure [T8]d. As with i{X), each evaluation here is proportional to the comoving number density 
of sources, ric, that give rise to the absorption times their average physical cross-section AA 
at that range of column densities. For A^hi < lO^^cm"^, the Ai{X) values rise steeply with 
decreasing A^hi- At A^hi ~ lO"'^'^'^ cm~^, however, the Ai{X) values flatten and are nearly 
constant until A^hi ~ lO^*^'^ cm~^. This simply reflects, of course, the 'inflection' in /{Nm, X) 
across the LLS regime. 

The flattenin g of /(A^hti X) as A^hj, decreases from y 10^° '^ has been recognized and dis- 
cussed preyioud^^jjO^Me^^ at higher redshifts, and also at 
2r ~ (ICorbelli &: Bandierall2002l ). This is expected to occur because of the transition in gas 
from a predominantly neutral state to a predominantly ionized plasma (jZheng fc Miralda-Escude 
2OO2I ) . This leads to a significantly lower cross-section for gas with A^hi < 10^° cm~^ than 
otherwise and hence an observed flattening in A£(X)/Alog A^hi- Detailed radiative trans- 
fer calc ulations of galaxies in high-reso l ution cosmological simulations have con firmed this 
picture (IFaucher-Giguere fc Kerea 1201 ll : iFumagalli et al.l I2OIII : lAltay et al.l 1201 ll ) . 



The steepening of f{N}ii,X) at A^hi ~ lO^^cm"^ with decreasing A^hi suggests a similar 
'phase transition'. Because A^hi = 10^^'^ cm~^ corresponds to Tg^ = 1, the obvious culprit 
for an inflection is the transition from the optically thin regime to optically thick gas. Qual- 
itatively, this may work as follows. For systems with A^hi lO^''''^ cm~^, the material is 
optically thin and one expects to receive roughly the same radiation field from the EUVB. 
In this case, the distribution of A ^ht values tha t result is dominated by the distribution of 
densities within the volume (e.g. ISchayd 1200 ll ). Our results indicate an A^hi distribution 
characterized by f{Niii,X) oc A^hj^. In the absence of radiation transfer effects (i.e. the 
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transition to optically thick gas), we may presume that the A^hi distribution would have 
shown a similar shape for A^hi > lO^^cm"^. Instead, as the gas becomes optically thick, the 
radiation field has been attenuated and one recovers a higher H I column density than in 
the optically thin limit. This non-linear effect naturally leads to the flattening of f{Ni{i,X) 
at A^Hi ~ lO^^cm"^. A full treatment of these effects within the context of cosmological 
simulations is clearly warranted. 



7.3. Evolution in /(A/^Hi, a:) 



Buoyed by the relative success of our analysis, we reco n sider the A^hi frequency distri- 
bution at 2 ~ 3.7, previously evaluated by iProchaska et al.l ( l2010l ). Table |9] lists the obser- 



vational constraints. For the Lyct forest (A^m = 10^^ — ip i^ s cm we adopt a constraint on 



the s lope of /(A^hi, X) based on line-fltting analysis (e.g. iKirkman fc Tytlerlll997l : iKim et al. 



2OO2I ). We allow for signiflcant uncertainty in this constraint since, as explored by these and 
other authors, issues such as continuum placement, completeness, and non- uniqueness of 
models to address line blending often arise, especially at higher redshifts. The normalization 
of f{N}ii,X) at these low column densities is set by the integrated o pacity r^g". At z ^ 4, 
there is signiflcant disagreement in the empirical estimates of r^g'^] Faucher-Giguere et al. 



( 2008bl ) rep ort a value Trff" ^ 0.8 while other au thors have published much lower values 
/'^Lyo ^ 0.7 lKim et al.lboOTl : iDall'Aglio et al.lboOS ). The difference in these values substan- 



'eff 



tially exceeds the estimated uncertainties, statistical and systematic. In our primary analysis 
we adopt the higher value but then comment on the implications of lower val ues. The other 



empi r ical constraints are estirn ates of f{Niii,X) from the SLLS and DLAs (jO'Meara et al. 



2OO7I : IProchaska fc Wolfd 120091 ) . and X^J^l and £(X) integral constraints (PWO09,POW10). 

Table [TO] summarizes the best-flt model for z = 3.7, error estimates, and the correlation 
matrix while Figure [21] presents a comparison of the best-flt model with the data. Similar to 
the z = 2A analysis, we flnd that a six parameter model (with identical A^hi pivots) provides 
a good description of the observational constraints at z = 3.7. More remarkably, the best-flt 
slopes all lie within la of the estimates at lower redshift. The only statistically signiflcant 
evolution lies in the normalization; ki2 is ~ 0.3 dex higher at higher redshift. 

The apparent lack of evolution in the shape of f{N}ii,X) between redshifts 2.4 and 3.7 
takes on additional interest when we consider the strong evolution in i{X) over the same 
timescale. As i{X) relates to the comoving number density of absorbers and their physical 
size, we can surmise from our results that although either (or both) of these two quantities 
must evolve with time, a change in the column density distribution of absorbers producing 
the LLS does not. This implies that if it is the change in physical size which is providing 
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the bulk of the evolution in ^(X), and if the LLS arise in galaxy halos, the distribution of 
the gas within the hal os is the same over cosmic time. The same conclusion is reached by 



Fumagalli et al.l ( 120111 ) in their simulations. 



7.4. Implications for the nature of LLS 



Finally, we consider these results in the context of current understanding on the origin 
of the LLS. At low redshift, it has been statistically established that syste nis with iVnT ^. 



10 cm are associate d with the circumgalactic medium (CGM) of galaxies (jProchaska et al. 



201 ll : iRudie et al.ll2012[ ). Specifically, galaxies of all luminosity show a very high covering 
fraction of significant Lya absorption to impac t parameters r ij a 300 k pc which may be 
defined as the extended CGM of these systems. iProchaska et al.l (120111 ) demonstrate that 
He for L > O.OIL* galaxies at 2; ~ is large enough that all of the strong Lya absorbers 
(> 0.3A) are associated with this extended CGM. At high z, one also observes a signific antly 
higher incidence of H I absorption in the CGM of L ^ L* galaxies (IRakic et al.ll201ll ) and 
a remarkably high incidence of optically thick gas in the CGM of galaxies hosting quasars 



(iHennawi fc Prochaskal 120071 ). One may infer, therefore, that gas with A^hi ^ 10 cm at 



z ~ 



2 primarily traces the extended CGM of high 2: galaxies. From numerical simu lations, 



one reaches similar conclusions (e.g. iFumagalli et al.l 1201 ll : Ivan de Voort et al.ll2012l ). 



Recent work in cosr nological siraulatio ns has also produced f{Ni{i,X) with numerous 
inflections. For example, lAltay et al.l (120111 ) provide a determination of f{Niii,X) at z = 3 
through a combination of line-fitting at low A^hi and projection at A^hi in the LLS and 
DLA regime. Their f{Niii,X) is consistent with those presented here at both z = 2.4 
and 2; = 3.7. Of particular interest is their observation that the presence of self shielding 
produces the flattening of /(Ahi,-^) in the LLS re gime, as discussed above. Unfortunately, 
their results are constrained to a single redshift. Ivan de Voort et al.l (120121 ) use the same 
set of simulations but focus instead on the gas fraction in halos, and the future history 
of LLS gas. They conclude that the majority of LLS gas is in cold-mode galaxy halo gas 
and will be incorporated in those galaxies' ISM within 1.2 Gyr. While they do not provide 
estimates of /(Ahi,-^) at other z, one may interpret their result along with our observed 
lack in evolution of /(Ahi,-^) to imply that to maintain the shape of /(A^hi,-^) at later 
times in the uni verse, cold gas must con tinue to be placed in galaxy haloes and be observed 
as LLS. Finally, IFumagalli et al.l (120111 ) also find that the majority of LLS are cold-mode 
gas. Notably, they follow their simulations through the redshift range 1.4 < z < 4 and find, 
in agreement with our results, that /(Ahi,-^) does not evolve with time. 
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8. Summciry and Future Work 

We have presented the first scientific results from our HST program to study Lyman 
limit absorption at redshift ;2 ~ 2. We have: 

• Determined i{z) and i{X) over the redshift range 1 < z < 2.6, finding a weighted 
mean of i{X)= 0.29 ± 0.05 for 2.0 < z < 2.5. 

• Stacked our WFC3 spectra to find a median (mean) value of the mean free path to 
lonizmg radiation of A^\^ = 243(252)/if2^ Mpc with an error on the mean value of 
±43/if2^ Mpc. 

• Shown either that the SED of 2; ~ 2 quasars is harder than those at 2; ~ 1 or that 
there are significant systematics in the SDSS selection function. 

• Constrained the evolution of A^\p with redshift, and estimating the breakthrough red- 
shift to be 2; = 1.6. The methodology behind this result is subject to the assumption 
that there are no significant inflections in the intrinsic QSO spectrum other than the 
one required to produce the observed Lyman alpha forest opacity. Including further 
inflections would produce very large uncertainties in the derived A^fp. 

• Shown that no 4-parameter (3 power-law) form can describe /(A^hi, ^) at 2; = 2.4, and 
provide a best-fit 5-parameter model. 

• Shown that f{Nm, X) does not vary significantly in shape between z = 2.4 and z = 3.7. 

While the current data sample are more than adequate to constrain the incidence of 
LLS and their contribution to /(iVHi,X) and A^\p at 2; ~ 2, they are of too poor resolution 
to make any statements on metal absorption in the LLS. Given their increased number over 
the DLA at these (and all) redshifts, the LLS likely contribute to a significant fraction of 
the metal budget. In the next paper in this series, we combine our HST results with those 
from Keck+ HIRES to address the metallicity of LLS at 2; ~ 2, a quantity which is currently 
poorly constrained. 
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Table 1. Quasar Sample and Continuum 



Object Name Plate" Fiber" Zem g 



Quasars Observed with WFC3/UVIS-G280 



J075547. 


.83+220450, 


.1 


1264 


163 


2, 


.319 


17, 


.58 


92.3 


0.4 


J075158, 


.65+424522, 


.9 


434 


555 


2, 


.453 


17, 


.77 


58.0 


0.6 


J084525, 


,84+072222, 


.3 


1298 


512 


2, 


.307 


17, 


,87 


49.8 


0.5 


J085045. 


.44+563618, 


.7 


448 


466 


2, 


.464 


18, 


,05 


65.7 


0.3 


J085316, 


.55+445616, 


.6 


897 


392 


2, 


.540 


18, 


.17 


74.8 


0.4 


J085417, 


.60+532735, 


.2 


449 


273 


2, 


.418 


17. 


.20 


113.4 


0.9 


J080620, 


,47+504124, 


.4 


1780 


181 


2, 


.457 


17, 


,57 


144.1 


0.0 


J083326, 


.82+081552, 


.0 


1759 


283 


2, 


.581 


17, 


.80 


90.9 


0.1 


J090938, 


.71+041525, 


.8 


1193 


296 


2, 


.444 


18, 


.01 


73.1 


0.2 


J094942, 


,34+052240, 


.3 


994 


245 


2, 


.282 


17, 


,75 


90.9 


0.5 


J105315, 


,89+400756, 


.4 


1435 


205 


2, 


.482 


18, 


,15 


49.8 


0.2 


J100541, 


.26+570544, 


.6 


558 


30 


2, 


.308 


18, 


.02 


90.5 


0.0 


J101120, 


.39+031244, 


.6 


574 


281 


2, 


.458 


17, 


.79 


127.7 


0.3 


J113550, 


,68+460705, 


.0 


1442 


1 


2, 


.496 


17, 


,67 


87.2 


0.4 


J114358, 


.52+052445, 


.0 


838 


591 


2, 


.561 


16, 


.95 


190.7 


0.3 


J110735, 


.58+642008, 


.7 


596 


385 


2, 


.316 


17, 


.38 


171.9 


0.1 


J111928, 


.38+130251, 


.0 


1605 


406 


2, 


.394 


18, 


.01 


56.9 


0.5 


J110411, 


,62+024655, 


.3 


508 


520 


2, 


.532 


18, 


,21 


41.5 


0.5 


J124831, 


.65+580928, 


.9 


1317 


577 


2, 


.599 


17, 


.74 


77.9 


0.2 


J125345, 


.49+051611, 


.3 


848 


438 


2, 


.398 


17, 


.99 


49.7 


0.4 


J125914, 


,85+672011, 


.8 


495 


494 


2, 


.443 


17, 


,93 


97.8 


0.2 


J121519, 


,42+424851, 


.0 


1450 


486 


2, 


.310 


17, 


,73 


124.0 


0.2 


J122015, 


.50+460802, 


.4 


1371 


327 


2, 


.446 


18, 


.01 


81.9 


0.4 


J122836, 


.05+510746, 


.2 


971 


448 


2, 


.450 


17, 


.88 


94.6 


0.0 


J123515, 


,84+630113, 


.4 


780 


546 


2, 


.383 


17, 


,38 


157.8 


0.4 


J133521, 


.96+454238, 


.2 


1463 


29 


2, 


.452 


18, 


.00 


114.3 


0.0 


J130055, 


.67+055620, 


.5 


849 


330 


2, 


.446 


17, 


.99 


95.6 


0.2 


J130240, 


.16+025457, 


.6 


524 


427 


2, 


.414 


17, 


.84 


117.5 


0.2 


J131123, 


,09+453159, 


.4 


1375 


459 


2, 


.403 


18, 


,20 


61.8 


0.0 


J131855, 


.75+531207, 


.2 


1041 


319 


2, 


.321 


18, 


.08 


70.3 


0.5 


J132312, 


.83+414933, 


.0 


1462 


482 


2, 


.440 


18, 


.06 


105.2 


0.4 


J132552, 


.17+663405, 


.7 


496 


114 


2, 


.511 


17, 


.92 


105.3 


0.3 


J133433, 


,88+035545, 


.2 


853 


220 


2, 


.583 


18, 


,14 


89.3 


0.2 


J133523, 


.82+463742, 


.1 


1463 


586 


2, 


.474 


17, 


.97 


86.7 


0.0 
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Table 1 — Continued 



Object Name 




Plate" 


Fiber" 


z 













J133646, 


,56+015741, 


,9 


528 


141 


2, 


,379 


17, 


,94 


112.0 


0.0 


J134211, 


,98+601525, 


.0 


786 


255 


2, 


.399 


17, 


,91 


101.2 


0.1 


J135412, 


,28+542100, 


,1 


1323 


213 


2, 


.294 


17, 


,98 


72.5 


0.3 


J135445, 


,66+002050, 


,3 


301 


385 


2, 


.504 


18, 


,07 


65.6 


0.0 


J135831, 


,78+050522, 


.8 


856 


469 


2, 


.455 


17, 


,17 


155.9 


0.4 


J141528, 


,47+370621, 


.2 


1643 


537 


2, 


.374 


17, 


,63 


96.3 


0.3 


J140028, 


,81+643030, 


,9 


498 


296 


2, 


.359 


18, 


,09 


56.4 


0.5 


J145453, 


,53+032456, 


Q 
.O 


Ooo 


iUo 


9 


.OOo 


1 7 




Qf^ A 
yD.4 


u.u 


J153335, 


.82+384301, 


.1 


1293 


30 


2, 


.529 


17, 


.50 


128.3 


0.2 


J153514, 


,65+483659, 


.7 


1167 


471 


2, 


.542 


17, 


.84 


108.0 


0.1 


J154042, 


,98+413816, 


,3 


1053 


306 


2, 


.516 


17, 


,46 


89.7 


0.0 


J162516, 


,42+294318, 


.3 


1421 


463 


2, 


.357 


17, 


,72 


70.6 


0.3 


J161003, 


,54+442353, 


.7 


814 


301 


2, 


.588 


18, 


,24 


46.7 


0.2 


J162548, 


,79+264658, 


.7 


1408 


281 


2, 


.518 


17, 


,34 


234.4 


0.0 


J165137, 


,52+400218, 


.9 


630 


560 


2 


.343 


17, 


,46 


109.1 


0.6 


J172409, 


,19+531405, 


.5 


359 


189 


2, 


.547 


17, 


,87 


52.8 


0.6 


J211157, 


,78+002457, 


.5 


986 


403 


2, 


.325 


17, 


,85 


94.0 


0.3 


J213629, 


,44+102952, 


,2 


731 


56 


2, 


.555 


17, 


,98 


59.1 


0.2 


J233823. 


,16+150445, 


.2 


747 


493 


2, 


.419 


17, 


,62 


108.9 


0.2 



Quasars Observed with ACS/PR200L 



J010657.94-085500.1 


659 


420 


2, 


.350 


18, 


,08 


2 


.6 


0.0 


J023359. 71+004938.5 


407 


600 


2, 


.522 


17, 


,85 


4, 


.0 


0.4 


J034024.57-051909.2 


462 


445 


2, 


.340 


17, 


,95 


3 


.9 


0.0 


J040241.42-064137.9 


464 


32 


2, 


.432 


17, 


,93 


3 


.0 


0.1 


J084619.53+365836.8 


865 


135 


2, 


.335 


18, 


,18 


3 


.3 


0.3 


J080413. 66+251633.8 


1205 


164 


2, 


.302 


17, 


,90 


1 


.9 


3.0 


J085952.36+350724.7 


1211 


506 


2, 


.373 


18, 


,22 


2 


.5 


0.6 


J092849.24+504930.5 


767 


41 


2, 


.349 


18, 


,14 


3 


.1 


0.1 


J102900.79+622342.0 


773 


314 


2, 


.470 


17, 


.78 


5, 


.9 


0.6 


J104321. 55+624127.3 


773 


177 


2, 


.355 


17, 


,74 


4, 


.3 


2.2 


J114311.61+651513.5 


597 


604 


2, 


.392 


17, 


,59 


4, 


.7 


0.0 


J121944.79+461015.3 


1451 


140 


2, 


.484 


18, 


,02 


2 


.5 


0.1 


J120729.75+042909.9 


842 


68 


2, 


.413 


18, 


,21 


1 


.5 


0.3 


J145554.30+521828.5 


1164 


288 


2, 


.348 


17, 


,77 


3 


.9 


0.7 


J165914.54+380900.7 


632 


120 


2, 


.343 


18, 


,12 


3 


.0 


1.2 
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Tabic 1 — Continued 



Object Name 


Plate" 


Fiber" 




9 






J161815.53+370103.9 


1056 


623 


2.496 


17.74 


4.1 


0.2 


J161538.88+392051.1 


1336 


137 


2.315 


18.11 


3.0 


0.0 


J230011.74-102144.4 


725 


295 


2.297 


18.26 


3.1 


0.0 



Note. — Quasars observed as part of our WFC3/UVIS-G280 and 
HST/ AGS survey for LLS absorption (Paper I). 

^SDSS-DR7 designations. 

''Absorbed continuum fitting parameters. The modeled continuum is of 
the form: /^°"*' = C/J'^'*'^'" (A/2500 A)«.We caution the reader that these 
models arc not meant to precisely describe the intrinsic spectral energy 
distributions of the quasars, but instead describes the flux attenuated by 
the Lyman series absorption (see the text). 



Table 2. LLS Models 



Quasar z^hs -^hi 



Quasars Observed with WFC3/UVIS-G280 



J075158.65+424522.9 


2.380 


16.50 




1.654 


17.10 




1.545 


17.35 


J075547.83+220450.1 


2.238 


16.10 




2.024 


17.25 




1.640 


16.80 


J080620.47+504124.4 


1.840 


16.95 




1.658 


17.15 




1.343 


16.80 




1.080 


17.05 


J083326.82+081552.0 


1.911 


16.40 


Quasars Observed with ACS/PR200L 


J010657.94-085500.1 


2.149 


16.40 




1.947 


18.10 


J023359.71+004938.5 


2.066 


16.60 




1.393 


16.35 


J034024.57-051909.2 


2.065 


18.15 


J040241.42-064137.9 


2.160 


17.10 




1.495 


16.85 




1.778 


16.65 



Note. — [The complete version of this table 
is in the electronic edition of the Journal. The 
printed edition contains only a sample.] 
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Table 3. WFC3/UVIS-G280 LLS STATISTICAL SURVEY 



Q 


uasar 




z 


em 


~r>2 
''start 


~r>l 
''start 


^r>0.5 
^start 


ZhhS 


log Nm 


JU ( oioo. 


1 /I 0/1 COO 
.DD+4z40zz, 


n 

.y 


o 


/ICQ 
.400 


1 


.576 


1 


.576 


1. 


.576 


1 


.576 


17. 


.60 


JU(£)04r 


QQ_LOon/i c^n 
.8o+zzU4oU, 


.1 


o 
z 


Q1 O 

.oiy 


1 


.400 


2, 


.026 


2. 


.026 


2 


.026 


17. 


.30 


JUoUOzU. 


/1 7_L^n/i 1 0/1 

.4 ( -|-0U4iz4, 


.4 




z, 


/I 

.40 ( 


1 


.400 


1 


.400 


1. 


.840 


1 


.840 


16. 


.95 


JUoOOZD. 


QO_Lr»ci ^P^o 
.oz+UolOOz, 


.u 




z. 


.Ooi 


1 


.400 


1, 


.400 


1. 


.400 












QA _Ln'70000 


Q 
.O 


o 

z. 


Qn7 
.oU ( 


2. 


.275 


2 


.275 


2. 


.275 


2 


.275 


18. 


.45 


JUo£)U40, 


A A _Lt=;^Q^1 Q 
.44+£)DODio, 


. / 


o 

z. 


.404 


1. 


.400 


1 


.400 


1. 


.400 










JUoOoiO, 


c^tt 1 /I /I c:^;i c 
.00+440010, 


.0 


o 

z, 


.04U 


1 


.903 


2, 


.117 


2. 


.117 


1 


.903 


17. 


.90 
























2 


.117 


17. 


.35 


JUo04i ( , 


.OU+Oc)z ( oO. 







z. 


/1 1 Q 
.4io 


1 


.400 


1 


.400 


1. 


.400 










JUyuyoo, 


71 1 ri/i 1 coc 

. ( i+U4i0Z0, 


Q 
.O 


z, 


AAA 
.444 


1, 


.400 


1 


.400 


1. 


.400 










Tno/io/io 


Q/i^n^oo/in 

.o4+U0zz4U. 


.6 




z, 


oco 

.zoz 


1, 


.400 


2 


.141 


2. 


.141 


2, 


.141 


17. 


.35 


JiUUo4i, 


.zO+0 ( U044, 


.0 


o 

z, 


.oUS 


2 


.239 


2 


.239 


2. 


.239 


2 


.239 


17. 


.75 


Ti ni 1 on 


QQ_l_nQ1 OAA 


.0 




z. 


.408 


1 


.634 


1 


.634 


1. 


.634 


1 


.634 


17. 


.60 


J iUOoiO. 


.oy+4UU / 00. 


.4 




z, 


A CO 
.4oz 


1 


.400 


2 


.050 


2. 


.050 


2, 


.050 


17. 


.35 


71 1 (7/1 1 1 

JiiU4ii, 


fio I no/i APiP^ 
.Oz+Uz4000. 


Q 


o 

z. 


CQO 

.Ooz 


1 


.613 


2, 


.062 


2. 


.495 


1 


.613 


18. 


.15 
























2. 


.062 


17. 


.40 
























2 


.495 


17. 


.15 


71 1 n7Qf^ 


.Oo+04zUUo. 


7 




z, 


Q1 A 
.oiD 


2, 


.093 


2 


.093 


2. 


.093 


2, 


.093 


18. 


.15 


7111 OOC 


.oo+ioUzOi. 


.U 


z, 


.oy4 


1 


.400 


2 


.149 


2. 


.149 


2 


.149 


17. 


.30 


71 1 Qc^cin 


.0O+40U (UO. 


.U 


o 

z, 


.4yo 


2 


.114 


2, 


.114 


2. 


.114 


2 


.114 


17. 


.55 


71 1 /I Q t^o 


.0z+U0z440. 




o 
z 


.001 


2, 


.139 


2 


.139 


2. 


.139 


2 


.139 


17. 


.50 


7101 CIO 

J izioiy. 


/lO 1 /10/1CC1 
.4z+4z4oOi. 


.U 


z, 


Q1 n 
.oiU 


1 


.400 


1 


.400 


1, 


.400 










71 ooni K 
J izzuio. 


.0U+4DUoUz. 


.4 




z, 


.440 


1 


.400 


1 


.400 


1. 


.891 


1 


.891 


16. 


.95 


71 OOQQA 


.U0+OiU(40. 




o 
z 


.40U 


1 


.400 


1 


.400 


1. 


.400 










J123515, 


.84+630113. 


.4 


2, 


.383 


2, 


.238 


2 


.238 


2. 


.238 


2, 


.238 


18. 


.50 


J124831, 


.65+580928. 


.9 


2, 


.599 


1, 


.400 


1 


.400 


2, 


.487 


2, 


.487 


16, 


.95 


J125345, 


.49+051611. 


.3 


2, 


.398 


2, 


.353 


2 


.353 


2. 


.353 


2, 


.353 


18. 


.55 


J125914, 


.85+672011. 


.8 


2, 


.443 


2 


.359 


2, 


.359 


2. 


.359 


2 


.359 


17. 


.95 


J130055, 


.67+055620. 


.5 


2 


.446 


2, 


.058 


2 


.058 


2. 


.058 


2, 


.058 


17. 


.60 


J130240, 


.16+025457. 


.6 


2 


.414 


1 


.917 


1, 


.917 


1. 


.917 


1 


.917 


18. 


.10 


J131123, 


.09+453159. 


.4 


2, 


.403 


1. 


.400 


1 


.400 


1. 


.809 


1. 


.809 


17. 


.05 


J131855, 


.75+531207. 


.2 


2, 


.321 


1 


.585 


1 


.585 


2. 


.080 


1. 


.585 


18. 


.00 
























2, 


.080 


17. 


.00 


J132312, 


.83+414933. 


.0 


2, 


.440 


2 


.267 


2 


.267 


2. 


.267 


2, 


.267 


18. 


.30 


J132552, 


.17+663405. 


.7 


2 


.511 


1 


.792 


2, 


.373 


2. 


.373 


1 


.792 


17. 


.60 
























2, 


.373 


17. 


.30 


J133433, 


.88+035545. 


.2 


2 


.583 


1 


.400 


2 


.269 


2. 


.269 


2, 


.269 


17. 


.45 
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Table 3 — Continued 



Q 


uasar 




z 


cm 


z^>2 

'^start 


'^start 


start 


2LLS 


log 


^HI 


J133521, 


.96+454238. 


,2 


2, 


.452 


2, 


.110 


2 


.110 


2, 


.110 


2, 


.110 


17, 


.65 


J133523, 


.82+463742. 


,1 


2, 


.474 


2, 


.327 


2 


.327 


2, 


.327 


2, 


.327 


17, 


.80 


J133646, 


.56+015741. 


,9 


2 


.379 


1 


.400 


1 


.400 


1, 


.400 










J134211, 


.98+601525. 


,0 


2, 


.399 


1 


.995 


1 


.995 


1, 


.995 


1 


.995 


17, 


.60 


J135412, 


.28+542100. 


,1 


2 


.294 


2 


.249 


2 


.249 


2, 


.249 


2 


.249 


17, 


.60 


J135445, 


.66+002050. 


,3 


2, 


.504 


1, 


.400 


1 


.400 


1, 


.400 










J135831, 


.78+050522. 


,8 


2, 


.455 


1, 


.895 


1 


.895 


1, 


.895 


1, 


.895 


18, 


.15 


J140028, 


.81+643030. 


,9 


2, 


.359 


2, 


.232 


2 


.232 


2, 


.232 


2, 


.232 


18, 


.45 


J141528, 


.47+370621. 


,2 


2 


.374 


2 


.124 


2 


.124 


2, 


.124 


2 


.124 


18, 


.45 


J145453, 


.53+032456. 


,8 


2 


.368 


1 


.861 


1 


.861 


1, 


.861 


1 


.861 


17, 


.50 


J153335, 


.82+384301. 


,1 


2, 


.529 


1 


.815 


1 


.815 


1, 


.815 


1, 


.815 


18, 


.35 


J153514, 


.65+483659. 


,7 


2, 


.542 


1 


.929 


1 


.929 


1, 


.929 


1 


.929 


18, 


.15 


J154042, 


.98+413816. 


,3 


2 


.516 


2 


.186 


2 


.186 


2, 


.186 


2 


.186 


18, 


.50 


J161003, 


.54+442353. 


,7 


2 


.588 


1 


.400 


1 


.400 


2, 


.331 


2 


.331 


17, 


.05 


J162516, 


.42+294318. 


,3 


2, 


.357 


1 


.400 


1 


.400 


1, 


.400 










J162548, 


.79+264658. 


,7 


2 


.518 


1 


.400 


1, 


.400 


1, 


.400 










J165137, 


.52+400218. 


,9 


2, 


.343 


1, 


.949 


1 


.949 


1, 


.949 


1, 


.949 


18, 


.10 


J172409, 


.19+531405. 


,5 


2 


.547 


2, 


.149 


2, 


.149 


2, 


.149 


2 


.149 


18, 


.00 


J211157, 


.78+002457. 


,5 


2, 


.325 


1 


.884 


1 


.884 


1, 


.884 


1 


.884 


17, 


.60 


J213629, 


.44+102952. 


,2 


2 


.555 


1 


.555 


1, 


.555 


2, 


.500 


1 


.555 


17, 


.60 
























2, 


.500 


17, 


.15 


J233823, 


.16+150445. 


,2 


2 


.419 


2 


.226 


2, 


.226 


2, 


.226 


2 


.226 


17, 


.75 



Table 4. ACS/PR200L LLS STATISTICAL SURVEY 



Q 


uasar 




z 


em 


yT>1 

^start 




U 


HI 


J010657, 


.94-085500, 


,1 


2 


.350 


1 


.947 


1 


.947 


18. 


,10 


J023359, 


.71+004938, 


,5 


2, 


.522 


1 


.200 










J034024, 


.57-051909, 


,2 


2, 


.340 


2, 


.065 


2, 


.065 


18, 


,15 


J040241, 


.42-064137. 


,9 


2, 


.432 


1 


.200 










J080413, 


.66+251633, 


.8 


2 


.302 


1 


.200 










J084619, 


.53+365836, 


,8 


2, 


.335 


1 


.200 










J085952, 


.36+350724, 


.7 


2, 


.373 


1 


.200 










J092849, 


.24+504930, 


.5 


2, 


.349 


1, 


.400 


1, 


.400 


17. 


,65 


J102900, 


.79+622342, 


.0 


2 


.470 


2, 


.279 


2 


.279 


17. 


,55 


J104321, 


.55+624127, 


,3 


2, 


.355 


1 


.200 










J114311, 


.61+651513, 


,5 


2 


.392 


2 


.101 


2 


.101 


18, 


,20 


J120729, 


.75+042909, 


.9 


2 


.413 


1 


.908 


1 


.908 


17. 


,70 


J121944, 


.79+461015, 


,3 


2 


.484 


2 


.171 


2 


.171 


18. 


,10 


J145554, 


.30+521828, 


,5 


2, 


.348 


1 


.200 










J161538, 


.88+392051, 


,1 


2 


.315 


2, 


.177 


2, 


.177 


17, 


,85 


J161815, 


.53+370103, 


.9 


2 


.496 


1 


.200 










J165914, 


.54+380900, 


,7 


2, 


.343 


1, 


.200 










J230011, 


.74-102144, 


.4 


2 


.297 


1 


.558 


1 


.558 


17. 


,85 
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Table 5. Incidence of LLS 



AX Az'' ml^g < z>' £{z) £{X) 



WFC3/UVIS-G280; Tg^L > 0.5 



[1.20,2.00] 
[2.00,2.60] 


22.9 
46.7 


8.0 13 1.85 
14.7 29 2.22 


1 62+°'*'* 

07+0.36 


0.57t^i^ 






WFC3/UVIS-G280; Tg^L 


> 1 




[1.20,2.00] 
[2.00,2.60] 


32.5 
52.5 


11.4 12 1.80 

16.5 25 2.19 


l.UO_Q 4Q 

-1 r 1+0.30 

^•^^-0.37 


0.37t°;J° 
0.48t°:?^ 






WFC3/UVIS-G280; Tg^^s 


> 2 




[1.20,2.00] 
[2.00,2.60] 


43.1 
56.2 


15.1 15 1.79 
17.7 17 2.20 


n QQ+O-25 

96+°-23 

U.»U_Q 2g 




^•'J'J-0.12 

■m+o-o'^ 

^•'^'^-0.09 


HST/ ACS; r^li > 2 


[1.20,2.00] 
[2.00,2.60] 


23.3 
17.4 


8.4 4 1.70 

5.5 5 2.16 


0.48t°i3 

01 +0-39 


Q -17+0.08 
^■-'^'-0.14 

29+°-i2 


Combined (with Ribaudo et al. 201 1~): r^h > 2 


[1.20,2.00] 
[2.00,2.60] 


108.0 
107.8 


39.0 29 1.77 
34.8 32 2.21 


0.69trJ 
0.92tri^ 


p. r,o + 0.06 



Note. — The cosmology assumed has ^Ia = 0.74, f2,„ — 0.26, 
and Ho = 72km s~^Mpc~^. For the combined values , we adopted 
the measurements from Table 7 of iRibaudo et all (|2011l ) for the z = 
[1.544, 1.947] and z = [1.947, 2.594] bins. We also converted their e{X) 
measurements to our cosmology with a simple scaling. 

^Total redshift survey path. 

^Number of LLS discovered in the survey path. Note that two LLS 
ocurring within w 10,000km of one another have been treated as 
one system. 
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'^Mean absorption redshift of the LLS. 
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Table 6. WFC3/UVIS-G280 STACKED SPECTRUM 



(A) 



601.92 0.46 0.09 

608.11 0.44 0.09 

614.30 0.42 0.09 

620.49 0.42 0.08 

626.68 0.43 0.08 

632.87 0.42 0.08 

639.06 0.42 0.08 

645.25 0.42 0.08 

651.44 0.42 0.08 

657.63 0.42 0.08 

663.82 0.42 0.07 

670.01 0.43 0.08 

676.20 0.43 0.07 

682.39 0.44 0.07 

688.58 0.44 0.07 

694.77 0.45 0.07 

700.96 0.46 0.07 

707.14 0.46 0.07 

713.33 0.44 0.07 

719.52 0.45 0.07 



'^Flux per A normal- 
ized to unity at Ar = 
1450A. 

''RMS (with 3(7 clip- 
ping) in the flux from 
a bootstrap analysis (see 
text). 

Note. — [The com- 
plete version of this table 
is in the electronic edi- 
tion of the Journal. The 
printed edition contains 
only a sample.] 
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Table 7. Stacked Spectrum Model Parameters and A^l 



Parameter 


Mean 


Median 


la c.l. 


2cr C.l. 


Ct 


1.03 


1.04 


0.99,1.05 


0.95,1.05 


SaT 


-0.37 


-0.33 


-0.60,-0.20 


-0.80,-0.20 


Lyman 


0.29 


0.28 


0.24,0.35 


0.21,0.44 


7r 


1.64 


1.73 


1.40, 1.80 


1.20,1.80 


K912(2stack) 


33.64 


33.93 


28.12,39.07 


24.42,47.15 


7« 


0.42 


0.00 


0.00,1.12 


0.00,2.55 


^912 


1.877 


1.890 


1.791,1.958 


1.713,2.018 


A^i2 (/,-iMpc) 


252.1 


243.4 


207.3,303.7 


177.3,353.1 
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Table 8. Revised A^\p Values from SDSS 



<Zq> 


\912 
'^mfp 

(Mpc-1) 


(Mpc-1) 


3.73 


54.9 


4.1 


3.78 


45.3 


3.6 


3.83 


45.4 


3.9 


3.89 


47.1 


4.2 


3.96 


38.9 


3.2 


4.08 


33.4 


3.0 


4.23 


27.8 


2.2 



Note. — Reevaluation of 
the A^jfp values using the 
stacked SDSS quasar spec 



tra from Prochaska et al 



(|2009f ). but update for cos- 
mology [Hq = 72km s~^ Mpc~^; 

= 0.26, = 0.74) and the 
defintion of the mean free path 
applied in this paper. 

''Based on a bootstrap analysis. 
There is additional systematic er- 
ror not included in this estimate. 



Table 9. /{N^lijX) Constraints 



Constraint 



log Nm 



Value" 



Comment 



Reference 



Constraints for f{N}ii,X) at z « 2.4 



Lya Forest 2.34 



SLLS 
DLA 



e{x)r>2 

^off 
\912 
^mfp 



2.51 



2.51 



2.23 
2.40 
2.44 



12 


.50- 


13, 


.00 


-11 


.19 


13 


.00- 


-13, 


.50 


-11 


.76 


13 


.50- 


-14, 


.00 


-12, 


.54 


14 


.00- 


-14, 


.50 


-13 


.30 


19 


.00- 


-19, 


.60 


-20, 


.63 


19 


.60- 


-20, 


.30 


-21 


.50 


20, 


.30- 


-20, 


.50 


-21 


.83 


20, 


.50- 


-20, 


.70 


-22, 


.29 


20, 


.70- 


-20, 


.90 


-22, 


.41 


20, 


.90- 


-21, 


.10 


-22, 


.88 


21 


.10- 


-21, 


.30 


-23 


.36 


21 


.30- 


-21, 


.50 


-23 


.60 


21 


.50- 


-21, 


.70 


-24, 


.33 


21 


.70- 


-21, 


.90 


-25, 


.00 


21 


.90- 


-22, 


.10 


-99.00t 


22, 


.10- 


-22, 


.30 


-99.00t 



> 17.49 
12.00-17.00 

12 - 22 



+0.05 
-0.04 
-f0.05 
-0.05 
+0.07 
-0.07 
+0.10 
-0.09 
+0.13 
-0.13 
+0.16 
-0.15 
+0.06 
-0.06 
+0.08 
-0.08 
+0.08 
-0.07 
+0.11 
-0.10 
+0.15 
-0.15 
+0.16 
-0.15 
+0.34 
-0.30 
+0.76 
-0.52 
-99.00 
-24.63 
-99.00 
-24.83 

0.30 ±0.07 
0.198 ±0.007 
243 ± 43 hj^ Mpc 



Recalculated for our Cosmology 



K02 



OPB07 



Only 30 systems total 
z = [2.3,2.7]; modest SDSS bias fsee Noterdaeme et al. 2009) ? PW09 



Converted to t^^" 



from Da- No LLS, no metals. 



This paper 
K05 

This paper 



Constraints for f{N}ii,X) at z ^ 3.7 



Lya Forest 
SLLS 



3.75 13.65-14.05 
3.58 19.00-19.60 
19.60-20.30 



-12.42^ 



-0.10 



-20.37l°:i2 

r,-, 99+0.14 



Read from Figure 4 
z = [3.1,4.5] 



KOI 
OPB07 



Table 9 — Continued 



Constraint 



log Nm 



Value'' 



Comment 



Reference 



DLA 



3.61 



20.30-20.50 
20.50-20.70 
20.70-20.90 
20.90-21.10 
21.10-21.30 
21.30-21.50 
21.50-21.70 
21.70-21.90 
21.90-22.10 
22.10-22.30 



-21.68tro^6 
-21.98tH? 

-22.45t°:}° 

-22.7i;«;JJ 

_OQ 00+0.17 

Z,O.ZO_Q j^g 



+0.18 
-0.17 



"23.53 
-23 93+0-24 

-24.53+°;f, 

-99.00+-^^° 
"99.00tl^^°° 



= [3.3,4.2]; modest SDSS bias (see Prochaska et al. 2009) PW09 





3 


70 


12.50-14.00 


-1.30 ±0.15 


Shallower than z=2, but very uncertain. 


See text 




3 


76 


> 17.49 


0.52 ±0.07 




POWIO 


^off 


3 


70 


12.00-19.00 


0.795 ±0.062 


Higher than other evaluations. 


FG08 


\912 
'^mfp 


3 


69 


12 - 22 


51 ± 5 /if2^ Mpc 




PWO09 



'^Effective redshift where the constraint was determined. 
^ f{Nm, X) constraints are given in log. 



References. — KT97: 



Kirkman fc Tvtlei (1997); KOI 



OPB07: lO'Meara et al.l (l20oJ; PW09: iProchaska et al 



Kim et all |200l[ ): K02: iKim et al.l (|2002[ ): K05: iKirkman etHI |2005l l: 
20091) 
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Table 10. Model and Correlation Matrix for the 5-Parameter /(A^hi, X) 



Param Best ct" 5ki2 5Pi2 



Model for /(iVni, X) a.t z !v 2.4 



kl2 


-9.52 


0.07 


1.00 


-0.91 


0.48 


-0.06 


0.01 


/3l2 


-1.67 


0.02 


-0.91 


1.00 


-0.70 


0.09 


-0.02 


017.5 


-1.07 


0.03 


0.48 


-0.70 


1.00 


-0.53 


0.11 


P20.3 


-1.71 


0.12 


-0.06 


0.09 


-0.53 


1.00 


-0.30 




-11.10 


7.54 


0.01 


-0.02 


0.11 


-0.30 


1.00 



'^Evaluated from the diagonal of the covariance matrix. 
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Fig. 1. — Comparisons between fit and input values for log A^hi (upper panel) and 2;lls (lower 
panel) for a mock sample of spectra designed to mimic the WFC3 data. The horizontal and 
vertical dashed lines in the upper panel correspond to r^^^ = 2. For most Tg^2 > 2 LLS, we 
can only obtain a lower limit on log A^hi and denote those systems with arrows. 
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Fig. 2. — WFC3/UVIS-G280 quasar spectra (black histogram) and the models (blue solid 
line) we have generated to assess Lyman limit absorption along the sightlines. The arrows 
mark the redshift and A'hi values for each modelled LLS. The dashed purple line traces our 
estimate of the continuum without LL absorption but includes Lyman series opacity. 
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Fig. 3. — Same as for Figure |2] but for the ACS/PR200L observations. Note that the x-axis 
has a logarithmic scahng. 




Wavelength (Aug) 
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g(z) 



N 



to 
b 



to 




Fig. 4. — These g{z) curves describe the survey size as a function of redshift. Specifically, 
g{z) expresses the number of quasars at a given z that were surveyed for LLS absorption to 
a limiting optical depth at = 9 12 A. They decrease with decreasing redshift because the 
presence of an LLS along a sightline often precludes the search for additional LLSs. The solid 
black curve shows the combined g{z) functions for the HST/ ACS and WFC3/UVIS-G280 
datasets for tqi2 > 2. 
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Fig. 5. — Incidence of LLSs per unit redshift i{z), as estimated from our iJST datasets. The 
WFC3/UVIS-G2803 data (black points) are shown for 3 hmiting optical depths riimit- The 
red (solid) points are for the ACS/PR200L sample and rnmit > 2. For each limiting optical 
depth, the data suggest only a modest (if any) evo l utioii in i{z) with redshift. The solid 



green line represents to the fit by ISongaila fc Cowid (j20ld ) for Tg{^ > 1. 
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(X) 




Fig. 6. — Incidence of LLSs per unit path length i{X) restricted to riimit > 2, as estimated 
from our combined WFC3/UVIS-G280 and HST/ ACS surveys. Sh own in sohd g r een is 
the i{X) curve (converted to the cosmology adopted in this paper) of iRibaudo et al.l (120111 ) 
derived from a single power-law fit to archival analysis of lower z HST data and high z obser- 
vations (jProchaska et al.ll2010l : ISongaila fc Cowid 120101 ). Their results are in good agreement 
with our estimates at 2; < 2.5 and together the data indicate mild evolution in i{X) at low 
z. 
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N 




Fig. 7. — The points (and error bars) show the incidence per unit absorption length i{X) 
for LLSs with hmiting optical depths rnmit = 0.5 — 2, for the WFC3/UVIS-G280 sample 
only. The red, blue, and green shaded regions express these same qua ntities. The orange 



shade d region shows an estimate for the incidence of DLAs taken from iProchaska &: Wolfe 
( 120091 ). Similarly, the yel low shaded reg i on gi ves the contribution of super Lyman limit 
systems (SLLSs) to i{X) Jo'Meara et allboOTh . Together, systems with Nm > 10^^ 
have i{X) ^ 0.18, implying a contribution of 50 — 80% to i{X)r>2- 



cm 
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f?^/fl450 




Fig. 8. — The black histogram shows the average quasar spectrum at 2; = 2.44 as estimated 
from a stack of the WFC3/UVIS-G280 spectra, each normahzed to unit flux at = 1450A. 
Note the strong decrease in flux at Ar < 900A due to the onset and cumulative effect of Lyman 
limit opacity. The gray shaded region shows the RMS at each pixel as measured from stan- 
d ard bootstrap analy sis. The solid red curve shows the radio-quiet average quasar spectrum 
of iTelfer et aD J2002h . smoothed to the resolution and dispersion of the WFC3/UVIS-G280 
stack. Their estimate of the intrinsic quasar flux, when normalized at Ar = 1450A, lies 
below the WFC3/UVIS-G280 stack at Ar ~ 1300A and would imply an effective opacity 
at Lya r^"" that is significantly smaller than previous estimates. We infer, therefore, that 
our WFC3/UVIS- G280 quasar cohort has an average SED that is bluer (harder) than the 
Telfer et al.l ( 12002| ) sample at far-UV wavelengths. 
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Fig. 9. — Upper panel: Allowed values for the normalization Ct and tilt Sot of the 



Telfer et al.l ( 120021 ) quasar SED (relative to the Te lfer et al. 2002 spec trum) that give the 



effective Lya opacity r^" measured aX z = 2.14 by IXytler et al.l (2004). The black asterisk 
denotes the model that minimizes from = 1080 — 1140 A; this model of the intrinsic 
spectrum is shown in the lower panel. We find that the WFC3/UVIS-G280 stack favors a 
bluer SED ((5q;t < 0) than the Telfer spectrum. Lower panel: The black histogram is the 
WFC3/UVIS-G280 stacked spectrum and the dotted curve is the Telfer spectrum scaled and 
tilted with the 'best' Ct and (5aT values (black asterisk in the upper panel). The solid red 
line shows this intrinsic SED attenuated by a constant r^"=0.19 for A < ALy^. The shaded 
region shows the range of models corresponding to the range of Ct and 5aT values in the 
upper panel. The gray dashed lines indicate the region where Da was evaluated. 
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Fig. 10. — Upper panel: Estimations of the effective Lyman series opacity t^™^" as a 
function of the tih in the Telfer quasar spectrum Sa^- Naturally, one estimates larger t^^^^^ 
values for steeper quasar SEDs. We estimate t^™'^" ~ 0.4 at z = 2.44. Lower panel: The 
black histogram shows the stacked WFC3/UVIS-G280 quasar spectrum with the error bar 
at Ar = 912A showing a 10% error estimate in the average observed flux at that wavelength. 
The solid red curves and the shaded region show the range of quasar SEDs that satisfy the 
constraint (Figure These are used in the estimates of t^™'''^ shown in the upper 
panel. 
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Fig. 11. — Upper panel: Estimated constraints on K9i2(-2stack) and 7^, our model parameters 
for the opacity of the universe to ionizing radiation at 2;stack = 2.44 (Equation [T2|) . These 
parameters are highly degenerate with one another. Lower panel: The black histogram is 
the WFC3/UVIS-G280 stacked quasar spectrum, normalized at Ar = 1450A. The purple 
dotted curve shows the best estimate of the intrinsic quasar continuum, i.e. the scaled and 
tilted Telfer spectrum as given in Table [71 The solid green curve gives the intrinsic SED 
attenuated by Lyman series opacity r^g^^^, constrained to match the stacked spectrum at 
912A to within 10% and assumed to decrease as (1 + z)"'^ . Lastly, the solid red curve is the 
complete model that also includes the effective Lyman limit opacity Tg^. 
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Fig. 12. — Each panel shows a crude estimate of the allowed Zgf^^ values corresponding to 
Tq^=1 from our analysis of the WFC3/UVIS-G280 stacked quasar spectrum, as a function of 
(upper) the quasar tilt 6aT and (lower) the assumed redshift evolution for t^™^°. The plus 
sign indicates values from the 'best fit' model shown in Figure [TTl Aside from Sa-r values 
near zero, we find 2912 = 1.85 — 2.1 implying A^\p~ 150 — 300/;.72^ Mpc. 
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Fig. 13. — Bootstrap analysis of A|^fp measured from the WFC3/UVIS-G280 stacked 
quasar spectrum (2;em=2.44) and using the techniques describe in § [5l The values peak 
at ^ 250hj2 Mpc with a non-Gaussian tail extending to larger values. 
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Fig. 14. — Estimates of the A^hi frequency distribution f{Niii,X) for systems with Nm > 
]^gi7.5 at 2 2.4. The measurements fo r the super Lyman hm it s ystems (SLLSs) and 
damped Lya systems (DLAs) were taken from lO'Meara et aLl (120071 ) and lProchaska fc Wolfe 
(I2OO9I ) respectively. The blue band is an estimate of f{Nm,X) for LLSs having Nm = 
_ iQi^cm"^ for an assumed power-law f{N}ii,X) oc A^^j^^ and constrained by the 
observed incidence of SLLSs and £(X)t->2 (this paper). We find /3lls = — 0.61q;3 (68% c.l.) 
for conservative estimates on the value of /(iVHi,X) at iVni = 10^^ cm~^ and allowing for 
the uncertainty in £{X)r>2- The dashed and dotted curves indicate the range of power-laws 
that satisfy the observations. 
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Fig. 15. — Upper panel: Tlie points witli error-bars sliow tlie measured offset Ai{X) 
in incidence of LLS i{X) from tlie Ti^^it > 2 estimate plotted at the limiting opti- 
cal depth. The curves show a series of model estimates for Ai{X) as a function of 
Tiimit for an assumed /(A^m = 10^^'^ cm~^, X) of 10~^^ cm^ (red), 10~^^'^cm^ (orange), 
and 10~^®cm^ (blue). We also assume a range of slopes /3pLLS for an assumed single 
power-law of —1 (solid), —1.5 (dotted), —2 (dashed) and —2.5 (dash-dot). Models with 
f{Nm = 10^^'^ cm~^, X) < 10~^^cm^ are ruled out by the observations. Similarly, mod- 
els with f{Niii = 10^^-^ cm"^, X) ^ 10^^*-^ cm^ prefer /3pLLS < -1-5. Lower panel: Per- 
mitted values for /3pLLS and /(Xhi = 10^^'^ cm~^, X) for an assumed single power-law 
(Equation [23]) that match the observations within the la uncertainty. Unless one adopts 
a large /(Xhi = lO^'''^ cm~^, X) value, the observations require a very steep power-law 
(/3pLLs < -1.5). 
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Fig. 16. — Estimates for the change in slope A^lls of the A^hi frequency distribution from 
Nhi ~ lO^^cm"^ to lO^^cm"^ as inferred from the incidence of LLS at 2; 2 in our HST 
survey. Results are plotted as a function of the incidence of r > 2 LLS, £{X)t->2- For our 
preferred value i{X)r>2 = 0.29 ±0.05, the data indicate a significant steepening of /(A^m, X) 
A/9lls > 1, as one transitions from optically thick gas to optically thin regions. 
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Fig. 17. — The curves trace the cumulative effective opacity for Lyman hmit absorption 
'^912 evaluated from ^912 = 1.85 to z^m = 2.44 (Equation [11]). Each curve assumes the 
identical /(A^hi,-^) distribution for A^hi > lO^^'^cm"^ as described by Figure [TH For A^hi < 
lO^^'^cm"^, the curves adopt a fixed power-law with /3 = —1.0 (dotted), —1.5 (dashed), 
—2.0 (dash-dot) and —2.5 (dash-dots). The red shaded region shows our best estimate for 
t}^i2 o'^^^ this redshift interval, as assessed from our A^£p analysis (§ I7.ip . This constraint 
strongly prefers a steep A^hi frequency distribution at these intermediate column densities, 
consistent with our analysis of the pLLS (Figure [T5l). 
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Fig. 18. — The A^hi frequency distribution function a.t z = 2.4, as modeled by a series of 
broken power-laws. The black curve shows the best-fit model, with parameters listed in 
Table [TOl This model was derived by minimizing the reduced calculated against the 
observational constraints shown in the Figure (see also Table |9]). This incl udes estimates 
of f(NuuX) for the Lya for est fKOl: Ixim et allbonil) . the SLLS fOPBOT IcTMeara et al. 



20071 ) ■ and the DLAs (PW09 IProchaska fc Wolfelboooh. Furthermo re, we adopted a model 



constraint on the power-law index in the Lya forest (jKim et al.ll2002[) and integral constraints 
on the opacity of the Lya forest (IKirkman et al.ll2005[ ). Lastly, we also adopted results on 
A^£p and i{X) estimated in this paper. We find that this simple (and assuredly non-physical) 
model for /(A^hi,-^) provides a good description of all these observations. 
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Fig. 19. — The data points trace the A^fp measurements as a function of redshift from this 
paper {z ~ 2.44; black) and the SDSS analysis of PWO09 {z ~ 4; green). The colored 
curves show the predicted evolution in A^fp for a series of assumed f{N}ii,X) models. Each 
has the functional form dervied for the IGM at z = 2.4 (the 5 parameter model given in 
Table [TOl) . We have allowed, however, the incidence of absorbers £{z) to evolve as (1 + z)'' 
with 7 = 1.5,2.0,2.5 for line colors brown, red, and blue respectively. Furthermore, each of 
the curves was forced to intersect the SDSS measurements ai z = 3.9. Our WFC3/UVIS- 
G280 measurement rules out 7 = 1.5 and favors 7 ~ 2 which coincides well with the redshift 
evolution in i{z) reported by iRibaudo et al.l (1201 ll ). The gray curve in the figure traces the 
Horizon of the universe with redshift. An extrapolation of the A^fp curves predicts that the 
'breakthrough' epoch occurs at z ~ 1.5 — 2. 
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Fig. 20. — Differential contribution of absorption systems to £{X)as a function of 
NB.i,di{X) / din Nm- Abrupt changes in this quantity mirror inflections in /(NjinX). The 
transitions away from the near constancy of this quantity across the LLS regime are likely 
linked to ionization properties of the gas, either with the onset of the gas being optically 
thick (A^Hi= 10^^-^cm"^), or the onset of self-shielding (A^hi= 10^°-^cm"^). 



- 82 - 



log f(NHi,X) 



to 



ON 



O 

OQ 



o 



to 




Fig. 21. — Similar to Figure [TSl but for z = 3.7 and, therefore, for a different set of obser- 
vationa l constraints (Table [9]). For r^", we adopted the estimate of iFaucher-Giguere et al. 
(l2008bl ) and the i{X) and A^\p measurements from POWIO and PWO09 respectively. The 
dotted gray line shows the best-fit model of f{N}ii,X) a.t z = 2.4. It has nearly identical 
shape but is offset to lower normalization by ~ 0.3 dex. 



